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A FOREWORD ON LAUNCH AND EXPERIMENT DATES 


When this Plan was written, the launch of SAGE was scheduled for January 25, 1979 
However, because of weather and spacecraft difficulties, the launch was delayed until February 
18, 1979 (1118 Eastern Standard Time) At that time the text of this Plan was already typeset, 
and therefore only minor changes have been made to the launch and correlative experiment 
dates stated in the body of this Plan 

At this time (February 20, 1979) the following guidelines on experiment dates can be 
given Correlative Experiment 2, originally scheduled for February 12-16 at White Sands, is 
now planned for March 11-12, with April 5-6 as backup dates Correlative Experiment 1, origi- 
nally scheduled for February 9 at Wallops Island, will also be rescheduled Although exact 
dates have not been set. Experiment 1 is now expected to occur after Experiment 2 (i e , after 
March 12) 

The dates of Correlative Experiments 3-6 are also subject to change, but specific plans 
have not yet been made 
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1. INTRODUCTION 


1.1 SAGE Measurements and Ground Truth Requirements 

The Stratospheric Aerosol and Gas Experiment (SAGE) is scheduled for launch in mid- 
February 1979 aboard the Applications Explorer Mission B (AEM-B) satellite of the National 
Aeronautics and Space Administration SAGE’s mission is to map vertical profiles of ozone, 
aerosol, nitrogen dioxide, and Rayleigh molecular extinction around the globe The ozone data 
are expected to extend from about 10-45 km, the aerosol data from cloud tops to about 35 km 
(plus occasional strong layers in the mesosphere), the nitrogen dioxide data from about 25-40 
km, and the Rayleigh molecular extinction data from about 15-40 km 

The SAGE instrument is a four-channel photometer that measures the intensity of sun- 
light (centered at wavelengths 0 385, 0 45, 0 60 and 1 0 ^.m) traversing the earth’s limb during 
spacecraft sunrise and sunset (See Appendix A for a description of the instrument ) In this 
manner it will measure vertical profiles of four-wavelength extinction, at the rate of about 30 
profiles per day Spatial coverage will extend from about 79° N to 79° S latitude (with some 
seasonal variation) and thus will complement the coverage (64° - 80° N and S) of the SAM-II 
stratospheric aerosol sensor on the Nimbus 7 satellite Figure 1 shows an example of latitude 
coverage versus time for a probable set of launch parameters 

SAGE’s four-channel extinction measurements will be numerically inverted to yield verti- 
cal profiles of ozone concentration, aerosol extinction (and inferred number density), nitrogen 
dioxide concentration, and total molecular density (When available, molecular density may be 
derived from the rawinsonde network and other sources, and then used as an input to the 
four-channel inversion process— if doing this improves the accuracy of the other derived param- 
eters ) The derived data will be archived and made available to the scientific community for use 
in a variety of studies However, before being released the data must be validated by compari- 
sons with correlative measurements made by other sensors of appropriate accuracy, resolution, 
and reliability To avoid confusion and the compromising of SAGE data integrity, the accuracy 
and resolution of correlative sensors should be well understood, if possible, they should be 
equal to or better than those expected from SAGE 

The anticipated performance of SAGE is indicated by Figure 2, which shows results of 
inversions using a typical SAGE inversion algorithm with simulated radiance data and simulated 
errors of the magnitude expected for a typical measurement cycle (Chu and McCormick, 1979) 
Although the accuracy and resolution achieved by SAGE are to a certain extent affected by lati- 
tude and by constituent concentrations, the models shown in Figure 2 give results that are 
representative of a broad range of conditions to be encountered by SAGE 

Note that the vertical resolution of each inverted profile is about 1 km This resolution is 
achieved by virtue of two factors (1) the radiometer’s narrow field of view, and (2) the 
sharply-peaked weighting functions for the limb-viewing geometry The limb weighting func- 
tions do in fact depend on constituent profiles, which can produce significantly poorer vertical 
resolution for certain conditions and heights Nevertheless, a vertical resolution of 1 km is a 
useful target specification for ground truth measurements 

Note also the error bars in Figure 2 These bars were derived by performing the inversion 
for ten different cases of simulated random errors and taking the standard deviation of the 
resulting set of solutions (The simulated random errors are based on SAGE expected 
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TANGENT LATITUDE 



FIGURE 1 SAGE LATITUDE COVERAGE AS A FUNCTION OF MISSION TIME 

Calculations assume an orbit inclination of 55°, a height of 600 km, and a launch time of. 1030 
EST, 25 January 1 979 







FIGURE 2 COMPARISON OF MODEL EXTINCTION PROFILES WITH RESULTS INVERTED FROM SIMULATED 
SAGE RADIANCE DATA 

Error bars show ± one standard deviation about the mean of inverted solutions for independent 
simulations with different random errors All simulations assumed probable errors of 3 arc seconds 
for pointing, and 0 5% (of signal) for radiometric measurement All radiometric data were quantized 
to 1 0-bit A/D accuracy 




performance, see Appendix A ) Hence, the error bars and the differences between model and 
inverted profiles denote the expected accuracy of SAGE-inferred constituent profiles and, there- 
fore, the desired accuracy for correlative measurements These desired accuracies are summar- 
ized in Table 1 


Table 1 

DESIRED ACCURACIES OF CORRELATIVE MEASUREMENTS 
(Based on SAGE Simulations) 


Ozone 

NO; 

Aerosol 

Molecular Density 

5%, 25—40 km 
10%, 10-25 km 
30%, 40—50 km 

30%, 25—40 km 
factor 3, 40—50 km 

10%, tropopause— 25 km 
30%, 25-35 km 
Plus noctilucent 
cloud occurrences 

5%, tropopause— 35 km 
10%, 35—45 km 


1.2 Overview of Ground Truth Plan 

A SAGE ground truth program has been developed on the basis of the above prere- 
quisites The program is designed to take advantage of established measurement capabilities 
and expertise at a number of worldwide locations At the same time it includes the develop- 
ment and use of sensors with needed measurement capabilities that cannot be provided by 
existing sensors 

A schedule of SAGE ground truth experiments to be supported by U S organizations is 
shown in Table 2 The sites for these experiments were selected because they permit a number 
of sensors for the SAGE constituents to be operated simultaneously from a single location or a 
tight cluster of sites Table 3 gives more information on the clustered sites that were con- 
sidered m developing the schedule of Table 2 Figure 3 shows the locations of the various clus- 
ters 


The intent of the schedule in Table 2 is to provide 

• A reasonably complete check of SAGE performance early in the mission 
(via the February-March 1979 experiments) 

• A check of any possible differences in SAGE sunrise and sunset perfor- 
mance (via the White Sands-Hollman experiment) 

• Ground truth support for a SAM-II/SAGE comparison (by means of the 
late spring or early summer 1979 Sondrestrom experiment) 
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Table 2 


SAGE GROUND TRUTH EXPERIMENTS SUPPORTED BY U S ORGANIZATIONS 
(To be supplemented by non-U S ground truth groups and others, as appropriate) 


Experiment 

1 

2 * 

3 

51 

5 

6 

Schedule 

(approximate) 

April 

1979 

March 11 
(sunset) 
March 12 
(sunrise) 
1979 

April 

1979 

Spring or 
Summer 
1979 

Fall 

1979 

Winter 

1979 

Cluster 

Cl 

C2 

C4 

C3 

Cl 

Cl 

Sites 

Wallops I 
Hampton 

Boulder 
White Sands 
Holloman AFB 

Fortaleza 

Natal 

Sondrestrom 
or Poker Flat 

Wallops I 
Hampton 

Wallops I 
Hampton 

Latitude 

Longitude 

37 - 38 N 
75 - 76 W 

33 - 40 N 
105 -107 W 

3 - 5 S 
35 - 38 W 

65 - 68 N 
51 or 148 W 

37 - 38 N 
75 - 76 W 

37 - 38 N 
75 - 76 W 

P-3 Lidar 

X 

X 

X 

X 

X 

X 

P-3 Spectrometer 

X 

X 

X 

X 

X 

X 

Langley 48" Lidar 
NCAR Lidar 

X 

X 



X 

X 

Dustsonde 


X 

X 

X 

X 

X 

O 3 Balloon 

X 

X 

X 

X 

X 

X 

O 3 Rocket 

K 

K(2),H(2) 

K 

K 

K 

K 

Datasonde 

X 

xxxx 

X 

X 

X 

X 

Murcray Balloon 
Interferometer 


X 





Pepin Balloon 
(photometer and 
spectrometer) 


X 


X 



UP* 

Noxon Spectrometer 


X 





Dobson/Canterbury 

X 

X 

X 

X 

X 

X 

NCAR Sabreliner§ 
Ames U-2 


X 


X 




"Includes correlative measurements for both a sunrise and a sunset SAGE scan 
^Combined SAM-II/SAGE correlative experiment 

*LIP (L1MS Instrument Package) is not formally a part of the SAGE Ground Truth Plan However, when LIP and SAGE 
measurements are sufficiently close in space and time, comparisons will be made A possible overlap may occur in Cold 
Lake, Alberta (55 N, 110 W) in February 1979 

{(Includes polar nephelometer, quartz crystal microbalance, Knollenberg optical particle counter, and possibly Dasibi O 3 sen- 
sor 

** 

Includes Ames aerosol impactor, Ames O 3 and NO sensors, Langley quartz crystal microbalance (planned), and possibly a 
Lazrus multifilter sampler 

K= Kruger optical rocket-borne O 3 sensor 
H= Hilsenrath chemiluminescent rocket-borne O 3 sensor 

NOTE Experiment dates are as scheduled on February 20, 1979 See also the Foreword on Launch and Experiment Dates 
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Table 3 


PRIME CANDIDATES FOR CLUSTERED SAGE GROUND TRUTH OBSERVATIONS 
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Table 3 (Continued) 


Cluster 

Number 


Instrument 


Investigator 


Ozone Balloon 

Perry, Bruton 

Ozone Rocket 

Perry, Bruton 

Datasonde 

Perry, Bruton 

Airborne Lidar* 

Fuller, 

McCormick 

Dustsonde 

Rosen, 

+ Oj Sensor 

Hofmann 


Agency or 
Institution 


Location Lat Long | Parameter 


Sondrestrom* 67 N* 

Sondrestrom 67 N* 

Sondrestrom* 67 N* 


5 S 35 



0 3 ,T 

35 W 

o 3 ,t 

35 W 

T,P,W 

35 W* 

A 


obson (Umkehr) 
zone Balloon 
erosol Lidar 


obson (Umkehr) 
zone Balloon 
erosol Lidar 
erosol Lidar 

erosol Lidar 


irborne Lidar 
uartz Microbalance 


ASA LaRC 
ASA LaRC 


5 S* 64 W* 
5 S* 64 W* 























Table 3 (Concluded) 


Cluster 

Number 

Instrument 

Investigator 

Agency or 
Institution 

Location 

Lat 

Long 

Parameter 

8 

Ozone Rocket 






mm 

8 

Ozone Balloon 







8 

Spectrometer 

Murcray 

Denver U 

Palmer* 

65 S* 

64 W* 


8 

Dustsondes 

Rosen, 

Hofmann 

U Wyoming 

McMurdo 

77 S 

165 E 

H 

9 

Ozone Balloon 

Kulkarm 

CSIRO 

Aspendale 

38 S 

145 E 

o 3 

9 

Dobson (Umehr) 

Kulkarni 

CSIRO 





9 

Aerosol Lidar 

Dilley 

CSIRO 

Aspendale 

38 S 

145 E 

A 

9 

Dustsonde 

Rosen, 

Hofmann 

U Wyoming 

Mildura 

34 S 

142 E 

A,T 

9 

Chemiluminescent 

Balloonsonde 

Galbally 

CSIRO 

Mildura 

34 S 

142 E 

NO, NOj0 3 

9 

Infrared 

Spectrometer 

Galbally 

CSIRO 

Mildura 

34 S 

14 1 E 

NO, NO^ O 3 


*Air-mobile 

^Ground-mobile 

^5/78, 10/78, 12/78, 2/79 Scheduled for LIMS verification (Dates listed in order of priority ) 

^On NCAR Sabreliner 

NOTE A = Aerosol, D = Density, NOj = Nitrogen dioxide, O 3 = Ozone, P = Pressure, T = Temperature, W = 
Winds 
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FIGURE 3 


LOCATIONS OF SOME POSSIBLE SAGE GROUND TRUTH CLUSTERS 

Sites 3a and 8a are also SAM II ground-truth sites (See, e g , Russell et al , 1 978 ) 



• Midlatitude, high-latitude, and low-latitude correlative measurements 

• Continuing checks of SAGE performance in different seasons 

Section 2 describes the sensors to be used in the United States-supported ground truth experi- 
ments Section 3 provides additional detail on scheduling and logistics, and Section 4 describes 
data-handling procedures 

Besides the U S -sponsored activities, both Japanese and European ad hoc ground truth 
groups have been formed, and efforts are being made to coordinate with Australian and Russian 
measurement teams These ad hoc groups include many established teams with a wide range of 
expertise in stratospheric constituent measurements The teams and their capabilities are 
described in Section 2 1 
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2. GROUND TRUTH SENSORS AND PROGRAMS 


2.1. SAGE Ad Hoc Ground Truth Groups 

A Japanese and a European SAGE ad hoc ground truth group have been formed The 
Japanese group is being coordinated by Professor Motokazu Hirono of Kyushu University in 
Fukuoka, the European group is being coordinated by Dr Robert Fantechi of the Commission 
of European Communities in Brussels Tables 4 and 5 show teams and sensors for each group 
There is also the possibility of establishing an Australian ground truth group, centered around 
the capabilities of CSIRO described in Table B-l of Appendix B, and possibly also a Soviet 
group under Professor K Ya Kondratyev of the Main Geophysical Observatory, Leningrad 

Each ad hoc group is now developing its own experiment schedule The role of the 
SAGE Experiment Team (SET)' in coordinating these activities will be limited, specifically, 
the SET will provide mission analysis data (i e SAGE measurement times and locations), 
specify desired data formats (see Section 4 2), and disseminate U S ground truth experiment 
schedules In addition to this coordination activity, however, the SET will participate fully in 
making use of Japanese, European, Australian, and Russian correlative data to validate the 
SAGE data 


2.2. Ozone Sensors 

2.2.1. Dobson Spectrophotometer 

A Dobson spectrophotometer 2 will be used to make ozone measurements at several 
ground truth sites (see Tables 2, 4, and 5) This instrument (Figure 4), by measuring selected 
ultraviolet wavelengths radiated by the sun, moon, or the zenith sky can be used to infer the 
quantity of ozone within the total air column above the instrument The result is expressed as 
a thickness of a layer of pure ozone at standard temperature and pressure 

Normally a Dobson station performs observations three times a day— in the midmorning, 
near local noon, and in the midafternoon At high-latitude stations fewer observations are 
called for during certain times of the year This will also be the case during any of the planned 
in situ ozone soundings from the various participating sites Because the accuracy of Dobson 
data deteriorates rapidly as the sun elevation angle decreases, twilight measurements are not 
practical 

Umkehr data (low-resolution vertical ozone profiles) will be obtained from certain Dob- 
son sites Efforts will be made to extrapolate data obtained at reasonable sun elevation angles 
to twilight points of tangency 


^SFT members are Or M P McCormick (leader) NASA Langley Research Center. Dr R A Craig Florida Slate University 
(deceased) Dr Derek M Cunnold Massachusetts Institute of Technology, Dr Gerald W Grams Georgia Institute of Technology Dr 
Benjamin M Herman University of Arizona Dr D t Miller British Meteorological Office Dr ») G Murcray University of Denver 
Dr T J Pepin University of Wyoming Dr Waller G Planet National Environmental Satellite Service and Dr Philip B Russell SRI 
International 

J U S Department of Commerce National Weather Service Ohv'ivtr s Manual Dobson Ozone Spei irophotonu ter Revised November 1 
1972 
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Table 4 


SAGE JAPANESE AD HOC GROUND TRUTH GROUP 


A Lidar Observation 


Team 

Objective 

Laser 

Site 

Department of Physics, 
Kyushu University 
(M Hirono) 

Aerosol Scattering 
10 - 40 km 

Ruby, Nd YAG 

* 

Fukuoka 

(33°37’N,130°26’E) 

Water Research Institute, 
Nagoya University 
(Y Iwasaka) 

Aerosol Scattering 
10 - 40 km 

Ruby, Nd YAG 

Nagoya^ 

(35°10’N,136°50’E) 

Department of Geophysics, 
Fohoku University 
(H Kamiyama) 

Aerosol Scattering 
10 - 40 km 

Ruby 

Zeo* 

(38°8’N,140°32’E) 


B Balloon and Spectrometer Observation (on Aircraft or Ground) 


Team 

Objective 

Instrument 

Site 

Upper Atmosphere Division 
Meteorological Research 
Institute 

1 no 2 , no, n 2 o, hno 3 
cfci 3 , cf 2 ci 2 , CC1 4 

1 Interference 
Spectrometer on 
Aircraft 

Tateno 

(36°3’N,140°8’E) 

(M Misaki) 

2 Aitken Particles, 
HCI, CH 4 

2 Balloon Sampling 


Atmospheric Research 
Institute 

Nagoya University 
(H Ishikawa) 

1 Aerosol Profile 

2 Ozone Profile 

1 Balloon Solar 
Occultation, 
Balloon Optical 
Particle Counter 

2 Balloon Solar 
Occultation 

Sanriku 

(39°8’N,141°49’E) 

Aerological Division 
Department of Observation 
lapanese Meteorological 
Agency 
(K Nyui) 

1 Ozone Profile 

1 Ozone Balloon 

1 Kagoshima 
(31°38’N,130°36’E) 
Tateno 

(36°3’N,140°8’E) 

' 

2 Ozone Column Content 

2 Dobson 
Spectrophotometer 

2 Tateno 
(36°3’N,140°8’E) 

Nugato University 

no 2 

Ground-Based 

Spectrometer 

Niigata 

(37°55’N,139°2’E) 


* 

Ruby observations made since October 1974, Nd-YAG observations expected to start in February 1979 
^Observations suspended at present, to be restored in 1979 
^Observations suspended at present, probably to be restored in 1979 
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Table 5 


SAGE EUROPEAN AD HOC GROUND TRUTH GROUP 


Team 

Objective 

Instrument 

Site 

ETH-Zurich 
(H Dutsch) 

1 Oj profiles, 
0-35 km 

1 O 3 sondes* 

1 Hohenpeissenberg 
(47°48’N,1 1°00’E) 
Payerne 

(46°49’N,6°57’E) 

Brussels 

(50 o 50’N,4 o 21’E) 
Elnas, Cagliari 
(39°13’N,9°08’E) 


2 O 3 profiles 
(higher but coarser) 

2 Dobson-Umkehr 

2 Arosa 

(46°47’N, 9°41’E) 


3 Total O 3 

3 Dobson 

3 Arosa 

(46°47’N, 9°41’E) 

Service d’Aeronomie, CNRS 
(M L Chanin, J Blamont, 
G Megie, P Aimedieu) 

1 03 ,Hi 0 profiles 
0-40 km, 
noctilucent clouds 

1 Dye lidars 

1 St Michel 
(44°N, 6 ° E) 
Verneres-le-Buisson 
(48°45’N,21°7’E) 
Heyss Island 
(80°5’N, E) 


2 O 3 , neutral density, 
aerosol profiles 
15-20 km (night) 
15-28 km (day) 

2 Long-life (—9 mos) 
variable-height balloon 
with chemiluminescent 
O 3 sensor, T,P sensors, 
plus U Wyo dustsonde 
or LRC aerosol sampler 
(possibly) 

2 Launch Praetoria 
(25°45’N,28°12’E) 
or Christchurch 
(43°33’S,172°40’E) 
drifts over wide area 


3 03 , 0 H profiles 

3 Rocket 

(Franco-Russian Program) 

3 Russia, India 


4 O 3 profiles, 
0-48 km 

4 Chemiluminescent 
balloonsonde 

4 Aire sur PAdour 
(43°42’N,0°15’W) 
and others 

Inst d’Aeronomie Spatiale 

Brussels 

(P Simon) 

1 O 3 profiles 
25-60 km 

1 Dropsonde solar 
integrating radiometer 
(Balloon [25-40 km] or 
rocket [ > 40 km] 
launch) 

1 UK Met Office rocket 
Possible rocket 
South Uist 
(57°N, 7°W) 

ONERA 

(J Laurent, A Girard) 

1 0 3 ,N0;,HAHN0 3 ,CFM 
profiles 

1 Solar IR grating spectrometer 
(balloon-borne) 

1 Aire sur PAdour 
(43°42’N, 0°15’W) 
Palestine 
(32°N,96°W) 

Inst Atmos Environ Res, 
Garmisch-Partenkirchen 

1 Aerosol profiles 
5-35 km 

1 Ruby lidar. 
Conductivity sonde 

1 Garmisch-Partenkirchen 
(47°30’N,1 1°05’E) 

(H Jaeger, W Carnuth, 
R Reiter) 

2 Total O 3 and 
profiles 

2 Matthews radiometer, 
O 3 sondes 

2 Garmisch-Partenkirchen 
(47°30’N,1 1°05’E) 
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Table 5 (Concluded) 


Team 

Objective 

Instrument 

Site 

AERE Harwell 
(A Eggleton) 

1 Oj horizontal 
variations near 
17 5 km max 

1 Dasibi UV absorption 
sensor on 
Concorde aircraft 

1 London-Washington DC 
(51°30’N, 0°10’W- 
38°55’N,77°00’W) 


2 Oj profiles 
5-17 5 km 

2 As above 

2 Points of ascent 
and descent 

CEC, Joint Research Center, 
Ispra 

(F Geiss) 

1 Aerosol profiles 

1 Ruby lidar 

1 Ispra 

(45°48’N, 8°36’E) 

Meteorological Office, UK 
(D Miller, J Gibbs, L Simmons, 
John Harries, NPL) 

1 Oj, H/) profiles 

1 1R emission sonde 
(J Harries, NPL) 

1 Bracknell 
(51°26’N, 0°46’W) 
or USA 


2 NOj, HNO 3 profiles 

2 SIBEX balloon package 
(J Harries, NPL) 

2 USA Sicily 


3 Aerosol profiles 
10-30 km 

3 Dye lidar 

605 nm, possibly 750 nm 

3 Bracknell 
(51°26’N, 0°46’W) 


4 0 3 

! 

4 Dobson spectrophotometer 

4 Bracknell 
(51°26’N, 0°46’W) 
Seychelles 
(4°36’S,55°30’W) 
St Helena 
(16°0’S,5°42’W) 
Lerwick, Scotland 
(60°09’N, 1°09’W) 


5 O 3 profiles 

5 Dropsonde solar radiometer 
(P Simon, IAS, Brussels) 
(UK Met Office rocket 
launch, see IAS, above) 

5 South Uist 
(57°N,7’W) 


6 Neutral density 

6 Mark III rawinsonde, 
Stratospheric Sounder 
on Tiros N for gridded 
analysis (with Finger, NOAA) 

6 Bracknell 
(51°26’N,0°46’W) 

Appleton Lab, SRC, 
Slough 

(L Thomas, T Gibson) 

1 Aerosol profiles 
5-35 km plus 
noctilucent clouds 

1 Dye lidar 

1 Winkfield 
(51'27’N, 0°43’W) 

Umv Koeln, FRG 
(A Ghazi, H Paetzold) 

1 O 3 H 3 O profiles 

1 Optical balloonsonde 

1 Koeln 

850°56’N, 6°56’E) 


2 Total O 3 

2 Dobson spectrophotometer 

2 Koeln 

850°56’N, 6°56’E) 


* 

Comparison of 5 different sondes at Hohenpeissenberg scheduled for April 1978 Dobson spectrophotometers will also be com- 
pared 


14 


























2.2.2. Canterbury Photometer 

The University of Canterbury (New Zealand) narrow- bandpass filter spectrometer, as the 
name implies, uses six narrow- bandpass interference filters to isolate wavelengths of interest 
These filters are sequentially rotated in front of the entrance aperture, thus allowing the pho- 
tomultiplier tube detector to measure intensities at each wavelength twice a second All six 
wavelengths are in the ultraviolet region of the spectrum and were chosen to match the stan- 
dard Dobson A, C and D wavelength pairs The radiation source can be the sun or the zenith 
sky The instrument design allows allowing the field of view to be changed from 2 3° , for use 
with the sun, to 4 6° for use with the zenith sky 

The Canterbury photometer, an extremely portable instrument, is being considered for 
acceptance as a standard method to measure total ozone concentrations A comparison of the 
Wallops Flight Center Dobson and Canterbury instruments has been in progress since October 
1977 The results thus far indicate that differences related to air mass and ambient temperature 
changes are not negligible, but that empirical corrections can potentially reduce the discrepan- 
cies to ± 2% of the Dobson values Note An extensive ozone photometer mtercomparison, 
involving the Dobson, Canterbury, Russian M-83, and the Canadian Brewer instruments, is 
being conducted An early effort to secure correlative support for SAGE from all these will be 
made 


2.2.3. Balloonborne ECC Ozonesonde 

The Electrochemical Concentration Cell (ECC) balloon-borne ozonesonde is a light- 
weight, compact, and relatively inexpensive instrument developed for measuring the vertical 
distribution of atmospheric ozone 3 An interior view of the sonde is shown in Figure 5 4 It is 
electronically coupled to a standard NOAA Radiosonde (Section 2 5) and thus also provides 
atmospheric pressure, temperature, and humidity in addition to ozone measurements The 
ECC ozonesonde is suspended approximately 22 meters below a 1200-gram balloon with the 
radiosonde suspended about one meter below the ozonesonde Figures 6 and 7 depict the flight 
configuration This combination, using the transmitter of the radiosonde, transmits data to the 
ground-based AN/GMD system 

The precision of the ECC ozonesonde is currently estimated to be within 10 to 12% ( 1 cr) 


2.2.4. Super Loki Optical Ozonesonde 

The Super Loki Optical Ozonesonde 5 is a rocket-launched payload, ejected at rocket apo- 
gee, that provides ozone profiles between 70 and 15 kilometers The ozone measurements are 
made during descent of the sensor, which is attached to a parachute The sensor provides an 
absolute measurement based on known values of the ozone absorption coefficients during day- 
light hours only Figure 8 shows the payload in the launch configuration with the standard 
Super Loki rocket used for the normal meteorological sounding systems The telemetry system 
operates on 1680 MHz, transmitting the data to a standard AN/GMD system equipped with a 
PCM adapter kit 

^Vehicles and Sensors of the UPN 607 Applications Sounding Rocket Program, October 1977, Preliminary 
4 

Instruction Manual, Electrochemical Concentration Cell Ozonesonde Mode! ECC-3A, January t 1977 

^Vehicles and Sensors of the UPN 607 Applications Sounding Rocket Program October 1977 Preliminary Instruction Manual Flec- 
trochcmical Concentration Cell Ozonesonde Cell Ozonesonde Model FCC-3a January I 1977 
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FIGURE 5 ECC OZONESONDE 



ECC sonde 


Sac* osonde- 
ozoncsonda 
connector cable 


Radiosonde 


FIGURE 6 ECC OZONESONDE HOOK-UP TO RADIOSONDE 
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Balloon 


high 

di ameter 


Parachute 


Ozonesonde 


Radiosonde 


Corner Reflector 


FIGURE 7 



4.57 m (15 ft.1 
2.44 m (8 ft.) 


Total length 
Is approximately 
25 m 



OZONESONDE BALLOON AND TRAIN 
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The sensor of the Super Loki Optical Ozonesonde is a four-channel filter-wheel UV pho- 
tometer A planar diffuser plate, located at the optical entrance of the sensor, intercepts incom- 
ing sunlight, while the varying solar angle of incidence is compensated for electronically The 
UV filters consist of interference filter elements to define the spectral bands and a common 
broadband UV filter that blocks unwanted radiation in the visible spectrum 

The precision of the optical ozonesonde is estimated at 5% or better between 25 and 50 
km This estimate is based on the properties of signal-to-noise ratios of the raw data for the 
ozone algorithm 6 The system accuracy is dependent on two factors (1) the error of the ozone 
absorption coefficient, and (2) the changes in the UV filter characteristics after calibration The 
accuracy is believed to be better than 10% The system can provide good data for solar zenith 
angles of about 80° or less, for angles exceeding 84°, the data become unusable 


2.2.5. Super Areas Chemiluminescent Ozonesonde 

The Super Areas chemiluminescent ozonesonde is a rocket-launched payload ejected at 
apogee, that measures the ozone distribution between 70 and 15 km as the sonde descends 
through the atmosphere on a parachute The ambient air is sampled by self-pumping, that is, a 
ballast tank, connected to the atmosphere by means of an inlet pipe, remains in pressure equili- 
brium with the increasing external pressure as the sonde descends The chemiluminescent 
detector and a photometer are oriented along the axis of the inlet pipe and enable continuous 
measurement of the ozone 

The chemiluminescent detector is an improved version of the one carried on balloon ozo- 
nesondes several years ago Its luminescence is proportional to ozone flux, which is the ozone 
concentration times the flow rate It can be shown that the flow rate is proportional to the pres- 
sure rise in the ballast tank, which is measured in flight Before flight a calibration is performed 
for each sonde This calibration simulates actual pressures, flow rates, and ozone concentra- 
tions expected during flight The measurement principle and early flight results have been 
described by Hilsenrath (1969, 1971) The present system flies on a Super Areas meteorologi- 
cal rocket, as shown in Figure 9, and utilizes the standard 1680-MHz AN/GMD telemetry sys- 
tem The sensor signals are pulse-code-modulated (PCM) and, therefore, require decommuta- 
tion for data processing 

An error analysis indicates precision of about 8% This is calculated from expected ran- 
dom errors due to "noise" and systematic errors or flight-to-flight errors Most of this error is 
related to the uncertainty entailed in establishing the sensor’s sensitivity from the calibration 
procedures Two flights conducted 13 minutes apart showed a repeatability of 6% This was 
derived from the average difference separating the two ozone profiles from a mean profile at 
one-kilometer intervals This result compares favorably with the computed 8% precision 

The absolute error in the experiment depends on the uncertainty of the ozone concentra- 
tion in the sonde calibration, undetected ozone losses, and nonlinearities in the measurement 
The uncertainty of the ozone concentration depends on the ozone monitor used in the calibra- 
tion This ozone monitor is calibrated at the National Bureau of Standards by measuring 
absorption at wavelength 253 7 nm The absolute measurement error is computed to be 12% 
(independent of the precision described above), which includes a 4% error associated with the 
ozone monitor used in the calibration Comparison of chemiluminescent rocketsonde 


^Design of Optical Ozonesonde for the Super Loki Dart Rocket Draft Report NASA Ooddard Space Flight Center Cireenbell Mary- 
land 
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FIGURE 9 SUPER ARCAS CHEMILUMINESCENT OZONESONDE 


measurements with balloon-borne and optical rocketsonde measurements shows differences no 
greater than 20%, usually on the order of 10%. Altitude resolution between 60 and 20 km is 
finer than 1/2 km. A comparison of the chemiluminescent and optical sondes with a sounding 
from LRIR on the Nimbus 6 satellite is shown in Figure 10. 


2.2.6. Ozone Lidar 

Recent work on the differential absorption lidar (DIAL) technique has led to the develop- 
ment of ozone lidar systems (e.g., Megie, et al., 1977). This technique uses a tunable laser to 
vary the output wavelength of the lidar so that it coincides with the wavelength of an ozone 
absorption line. The absorption of the laser radiation by ozone molecules thereby provides a 
method for determining the vertical profile of ozone molecules. Lidar echoes are recorded as a 
function of range for a wavelength corresponding to the center of an ozone absorption line, as 
well as for a nearby wavelength that is not absorbed by ozone or other atmospheric molecules. 
The system described by Megie et al. obtained ozone profiles from the 18- to 28-km altitude 
interval with approximately 1.2-km altitude resolution-using ultraviolet wavelengths of 308 nm 
and 303.7 nm generated by a frequency-doubled rhodamine-6G dye laser. 
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COMPARISON - 29 JULY 1975 
WALLOPS IS VA 



MIXING RATIO {^gm/gm) 


FIGURE 10 COMPARISON OF OZONE MEASUREMENTS MADE BY CHEMILUMINESCENT 

ROCKETSONDE, OPTICAL ROCKETSONDE, LRIR SATELLITE SENSOR, AND ECC 
BALLOONSONDE 
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2.2.7. PAM II 


The PAM II (Preliminary Aerosol Experiment) experiment constructed by Dr Pepin is 
planned for launch on the Air Force P78-1 satellite on February 19, 1979 PAM II is a solar 
photometer that operates like the SAGE and SAM II experiments in that it measures solar 
extinction at spacecraft sunrise and sunset events However, instead of scanning the solar disk, 
the PAM II instrument points to the radiometric center of the disk and measures the total sig- 
nal from the full disk-in a manner simlar to the measurement that was made using SAM on 
the ASTP flight (Pepin 1977) 

PAM II has three optical channels at wavelengths of 0 43, 0 60 , and 1 0 /im Its mission 
is to map vertical profiles of aerosols and ozone by scanning the atmosphere during spacecraft 
sunrise and sunset events Because P78-1 will fly in a polar, sun-synchronous high-noon orbit, 
the latitude bands covered by the PAM II experiment will be restricted to 63 to 82 0 N and S 
latitude 

Because SAGE’s latitude coverage is expected to extend from the equator to 79° N and S, 
the coverage of the PAM II and SAGE sensors will overlap partially At the times overlap 
occurs there will be opportunities to compare measurements made by the two independent 
remote-sounding systems The PAM II ozone measurements will be available for extending the 
study of the SAGE vertical ozone profiles to higher latitudes 


2.2.8. Other Ozone Sensors 

A variety of sensors is available for in situ ozone observations on aircraft platforms It is 
likely that most measurements of this type will be made with the Dasibi ultraviolet absorption 
instrument, which monitors ozone on a continuous basis by means of a long-pass UV absorp- 
tion cell The light source in this system is a 253 7-nm Hg lamp To obtain corrections for 
lamp or electronic drift, the instrument continuously compares the signal from the sample cell 
with that of an identical reference chamber The gas stream flowing through the reference cell 
is first scrubbed of all ozone The difference between the signal from the sample chamber and 
its counterpart from the reference cell then constitutes a measure of the absolute ozone con- 
centration 


2.3. Aerosol Sensors 
2.3.1. Dustsonde 

Figure 11 shows a schematic drawing of the University of Wyoming balloon-borne dust- 
sonde that is planned for ground truth in the SAGE program Its mode of operation is as fol- 
lows Air sampled during balloon ascent and parachute descent is pumped at approximately 
0 75 1/min in a well-defined stream through the focal point of the condenser lens in the 2 5-hter 
scattering chamber, where the individual stratospheric aerosol particles scatter light into the 
microscopes The light pulses that can be observed with the microscope are detected and 
amplified by the photomultipliers By means of pulse height discrimination and careful labora- 
tory calibration with aerosols of known size and index of refraction, the integral concentration 
of aerosol particles with radii greater than 0 15 and 0 25 /im can be determined We will refer 
to these integral concentrations as N 15 and N 25, respectively 
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Two photomultipliers are used to enhance the signal-to-noise ratio by counting only coin- 
cident events from the two detectors The background noise for the system is mainly due to 
Rayleigh scattering from air molecules in the chamber at low altitude and from cosmic ray scin- 
tillation in the photomultiplier glass at high altitude The requirement for coincidence of events 
from the two detectors removes the contribution due to the cosmic ray scintillation The back- 
ground is measured approximately every fifteen minutes during the flight by having filtered air 
pass through the chamber The background produced by the Rayleigh scattering is negligible 
above a 10-km altitude Below this altitude the measured corrections for the background are 
employed The dustsonde is also equipped with rawinsonde temperature elements for recording 
the vertical temperature profile 

Resolution and Accuracy 

Figure 12 shows a typical dustsonde-measured profile of particle number density Note 
that below about 28 km the vertical resolution is better than 1 km in both particle size chan- 
nels Professor Rosen has performed an analysis of the accuracy of dustsonde measurements 
(eg Hofmann et al , 1975) The major sources of error are counting statistics and possible 
variations in the refractive index of particles The counting method and sensor channels are 
designed to minimize these errors, for stratospheric heights below 25 km they result in typical 
uncertainties of about 8% for both Channels I and II (r ^ 0 15 and 0 25 respectively) 
Above 25 km these errors tend to increase significantly because of poorer counting statistics 
and less accurate measurement of the sampled air volume 

Dustsonde measurements can be converted to an estimated particulate 1 0-jim extinction 
coefficient by using an assumed refractive index and a two-parameter size distribution fitted to 
the two-channel dustsonde data on N 15 and N 25 Figure 13 shows the dependence of the 
conversion ratio on optical model properties 

[The size distribution functions and refractive indices shown have been derived from 
measurements by various investigators-e g Hofmann et al (1975), Toon and Pollack (1976), 
Harris and Rosen (1976), Swissler and Harris (1976) However, to generate the complete range 
of values shown for each curve, parameters were varied, sometimes beyond the range of obser- 
vations Note that observations of N 15 /N ^ less than 2 are very rare, and average values for 
stratospheric layers several km thick are typically between 3 and 5 ] 

In a given dustsonde measurement the channel ratio, N 15 /N i5 , is known, but the particle 
size model and refractive index can in general only be estimated on the basis of previous meas- 
urements Thus, the conversion ratio uncertainty is given by the vertical spread in the curves 
above the measured value of N 15 /N ^ Figure 14 shows the one-standard-deviation spread for 
the different aerosol compositions As can be seen, the uncertainty in converting a two-channel 
dustsonde measurement to 1 0-pim extinction is thus about ±25% if particle composition is 
unknown, and about ±15% if the refractive index is known to be one of the two values shown 
in Figure 13 (1 e, either silicate or aqueous sulfuric- acid composition) A similar conclusion 
was obtained by Pepin and Cerm (1977) 


2.3.2. Airborne Lidar 

An airborne lidar for SAGE and SAM-II ground truth measurements has been developed 
under the direction of William Fuller at NASA Langley Research Center (LRC) The lidar 
design is based on a study by Evans (1977), with appropriate modifications by LRC personnel 
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FIGURE 1 2 EXAMPLE OF TWO-CHANNEL DUSTSONDE MEASUREMENT OF PARTICLE NUMBER DENSITY, 
OBTAINED AT MCMURDO STATION (78° S, 167° E) DURING NONVOLCANIC CONDITIONS (12 
JANUARY 1973) 


(a) Number of particles with radius >0 15 /urn (N 15 , measured by channel I) Smooth curves are 
lines of constant mixing ratio (particle number per mg of air) (b) Ratio of numbers of particles 
measured by channels I and II N 2 s is number of particles with radius > 0 25 /urn. and is measured 
by Channel II Ratio data are plotted at heights where Channel II data are available N 15 data are 
interpolated at these heights in computing ratios 
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Table 6 shows the design parameters, Figure 15 a photo of the lidar itself (Although it 
currently uses only a ruby laser, plans call for the Nd-YAG laser and associated detector to be 
added in time for at least some of the SAGE correlative measurements) Figure 16 illustrates 
the organization of the airborne lidar project 

The platform for the airborne lidar is the P-3 aircraft of the NASA Wallops Flight Center 
Resolution and Accuracy 

Measurements with the NASA Langley 48" ruby lidar, the SRI 16" ruby-dye lidar, and the 
NCAR ruby lidar, among others, have shown that stratospheric aerosol measurements can be 
made with a vertical resolution of 1 km or better (up to about 30 km) by accumulating photons 


Table 6 

DESIGN PARAMETERS OF AIRBORNE LIDAR 


Transmitter 

Ruby 

Nd YAG 

Wavelength (/xm) 

0 6943 

1 06 

Energy per Pulse (J) 

1 0 

05 

Repetition Rate (pps) 

1 0 

20 

Pulse Width (n sec) 

30 

20 

Beam Divergence (mr) 

1 0 

1 0 

Beam Diameter (cm) 

8 

76 

Receiver 

Diameter (cm) 

36 

36 

Field of View (mr) 

2 

2 

Filter Bandwidth (A) 

10 

10 

Optical Eff to PMT 

0 35 

0 35 

PMT Quant Eff 

0 10 

0 03 

Skylight Background* [w/(m 2 srA)] 

2xl0 -4 

1 3xl0~ 5 

Data Acquisition 

Bandwidth 

1-2 5 mHz 


ADC Rate 

10 mHz (max) 


ADC Resolution 

10 Bit 


ADC Memory 

2048 Words 


Computer Memory 

32 K 16-Bit Words 

Magnetic Tape 

45 IPS, 800 CPI, 9 Track 


*For zenith-viewing lidar flying above 6 km, with sun near horizon 
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FIGURE 1 5 THE NASA LANGLEY AIRBORNE LIDAR SYSTEM 




FID - SIB IRD - AMRB IRD - AMRB IRD - DDAS IRD - AMRB 


FID - FLIGHT INSTRUMENTATION DIVISION 
IRD - INSTRUMENT RESEARCH DIVISION 
AMRB - AEROSOL MEASUREMENTS RESEARCH BRANCH 


DDAS - DIGITAL DATA ACQUISITION SECTION 
SIB - SPACECRAFT INSTRUMENTATION BRANCH 


FIGURE 16 ORGANIZATION OF AIRBORNE LIDAR PROJECT 








for reasonable amounts of time The accuracy of the particulate backscattering coefficients 
derived from such measurements is a strong and complicated function of the laser wavelength, 
other Ildar parameters, skylight background, aerosol concentration, the proximity of the nearest 
radiosonde sounding, the validity of normalization procedures, and even the uncertainty in the 
ozone vertical profile To evaluate this accuracy for realistic situations we have developed a 
computer program that simulates the measurement and data analysis process, as shown in Fig- 
ure 17 (A parameter shown in Figure 17 is the scattering ratio, R, a central quantity derived 
in the analysis of stratospheric lidar data It is defined as R = (B p + B g )/B g , where B p and B g 
are respectively the particulate and gaseous backscattering coefficients ) 

At each appropriate step of the simulation the program computes the relative uncertainty 
in each derived quantity by using an analytical expression The sources of error include (1) 
signal measurement error, (2) molecular density uncertainty, (3) aerosol and ozone transmis- 
sion uncertainty, and (4) normalization uncertainty As a check on the analytical expressions 
for error propagation, random number generators (symbolized by circles in Figure 17) are used 
to inject random errors from sources (l)-(3) at appropriate points of the simulation [Error 
(4), normalization, affects the entire derived profile in a systematic way and should not be 
simulated by different random errors at each data point ] Table 7 lists the sizes of the error 
sources used in the simulations (Justification for the chosen error sizes is given by Russell et 
al , 1976a,b ) 

Figures 18-20 show the results of using the program to simulate airborne lidar measure- 
ments for different latitudes and aerosol conditions 


2.3.3. Ground-Based Lidar 

2.3. 3.1. NASA Langley 48" Ruby Lidar System 

Shown in Figure 21 is the Langley Research Center’s 48" laser radar (lidar) rystem It 
comprises two temperature-controlled lasers (ruby and neodymium-doped glass) mounted on 
either side of an f/10 Cassegrain telescope consisting of a 48-inch-diameter f/2 all-metal pri- 
mary and a 10-inch diameter secondary mirror A schematic of the system is displayed in Fig- 
ure 22 The detector package output is recorded by a high-speed data acquisition system Ana- 
log signals are amplified and bandwidth-limited, digitized at a 10-MHz rate with 8-bit accuracy, 
and then recorded on magnetic tape Pulse count data are amplified, discriminated, counted at 
a 200-MHz rate, and also stored on magnetic tape Altitude resolution is obtained by using the 
variable (1-, 5-, or 10-microsecond) bin widths that are available A 16K-word-storage com- 
puter is used to control the data acquisition system and provide data processing An X-band 
microwave radar, boresighted with the laser system axis, is used to ensure safe operation in the 
atmosphere A rotating shutter reduces laser fluorescence after Q-switching The entire system 
is mobile and can scan in elevation and azimuth at a slew rate of 1° per second 

Presently the 48" system is being updated for simultaneous two-channel measurements 
and the addition of new ten-bit accuracy analog-to-digital converters In addition, a real-time 
graphics display system is being added 
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FIGURE 1 7 SIMULATION PROCEDURE FOR EVALUATING LIDAR MEASUREMENT ERRORS 

X L is the lidar wavelength Each simulation is also performed for the SAM ll/SAGE 
wavelength, 1 0 /xm Circles symbolize random number generators that inject 
simulated errors into derived quantities at appropriate steps of the computation 
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Table 7 


ASSUMED SIZES OF ERROR SOURCES IN LIDAR DATA ANALYSIS 



Notes 

B, I = Detector output resulting from background light and internal noise, respectively 

T 3 ,r p = One way optical thickness of ozone and aerosol particles, respectively, between nor- 
malization altitude and altitude of analysis 

Rmdx = Maximum scattering ratio in Ildar profile being analyzed 

X L = Lidar wavelength 

b = Exponent of power-law approximation to wavelength dependence of particulate back- 
scattering between 0 69 /^m and A[_ F° r most practical purposes b = I 8 

* 

Assumes radiosonde density profile available within about 100 km and 6 hours of lidar 
measurement, and no intervening frontal activity 

^Assumes ±20% uncertainty in T} and ±50% uncertainty in r p 

*Based on typical aerosol concentration present at height of minimum mixing ratio in a long 
series of nonvolcamc and postvolcanic duslsonde measurements at Laramie, Wyoming, and 
on several ruby lidar/dustsonde comparison experiments The error sizes shown apply to 
the case in which the scattering ratio profile is normalized to force its minimum to equal 
the value expected from previous dustsonde measurements, rather than the value 1 00 
commonly used If instead R m , n = 1 00 is forced, the expected errors become asymmetric 
(always negative) and roughly twice as large 
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ALTITUDE — km ALTITUDE 



0 40 80 0 40 80 


RELATIVE UNCERTAINTY IN PARTICULATE BACKSCATTERING — percent 

FIGURE 1 8 LIDAR MEASUREMENT UNCERTAINTIES AND SIMULATED MEASUREMENTS FOR 
THE W-48 (LARAMIE, WYOMING) MODEL ATMOSPHERE 

(a) Model scattering ratios, simulated measurements (dots), and expected error 
bars (b) Model backscattering coefficients, simulated measurements (dots), and 
expected error bars (c),(d) Relative uncertainty in particulate backscattering, 
showing contributions by source 
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ALTITUDE — km ALTITUDE 


10 12 14 (Ruby) 




COEFFICIENT — (m-sr)' 1 




RELATIVE UNCERTAINTY IN PARTICULATE BACKSCATTERI NG — percent 

FIGURE 19 LIDAR MEASUREMENT UNCERTAINTIES AND SIMULATED MEASUREMENTS FOR 
THE P-10 (ALBROOK, PANAMA) MODEL ATMOSPHERE 

(a) Model scattering ratios, simulated measurements (dots), and expected error 
bars (b) Model backscattenng coefficients, simulated measurements (dots), and 
expected error bars (c),(d) Relative uncertainty in particulate backscattenng, show- 
ing contributions by source 
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ALTITUDE — km ALTITUDE 


1 00 1 10 1 20 (Ruby) 




COEFFICIENT — (m-sr)" 1 




RELATIVE UNCERTAINTY IN PARTICULATE BACKSCATTERING percent 

FIGURE 20 LIDAR MEASUREMENT UNCERTAINTIES AND SIMULATED MEASUREMENTS FOR 
THE L-3 (LONGREACH, AUSTRALIA) MODEL ATMOSPHERE 

(a) Model scattering ratios, simulated measurements (dots), and expected error 
bars (b) Model backscattering coefficients, simulated measuremens (dots), and 
expected error bars (c),(d) Relative uncertainty in particulate backscattering, show- 
ing contributions by source 
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THE 48 INCH (1 .22-METER) NASA LANGLEY RESEARCH CENTER LIDAR SYSTEM 
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FIGURE 22 SCHEMATIC DIAGRAM OF THE 48-INCH (1 .22-METER) NASA LANGLEY 
RESEARCH CENTER LIDAR SYSTEM 
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2.3. 3. 2. NCAR 60" Ruby Lidar System 


The National Center for Atmospheric Research (NCAR) in Boulder, Colorado has a 
ground-based ruby lidar that has been used for stratospheric aerosol measurements at various 
times since about 1969 The system is built around a 60-inch-diameter searchlight mirror It 
has a pulse energy of 1-2 J and a pulse duration of 60 ns The system’s field of view is fixed at 
the vertical 

This system will be operated by Dr Frederick G Fernald of Denver University and Mr 
Charles Frush of NCAR as part of the SAGE ground truth program Specifically, it was used in 
the Laramie-Boulder practice comparative experiment in September 1978 and will be used in 
the White Sands-Boulder correlative experiment in February or March 1979 (See also Table 2 
and Sections 3 1 and 3 2 4) 


2.3. 3. 3. CNRS Heyss Island Lidar System 

Dr M L Chanin of France’s Centre National de la Recherche Scientifique (CNRS) will 
conduct a program of noctilucent cloud observations by dye lidar techniques from Heyss Island 
(80 5°N), where a lidar station has been established since 1975 for sodium measurements The 
observations, part of a joint Franco-Soviet program, are planned to start prior to the SAGE 
launch Dr Chanin has agreed to make available to the SAGE and SAM II ground truth pro- 
grams any NLC observations that are near SAM-II or SAGE scans in both space and time 


2. 3. 3. 4. NOAA Point Barrow Dye Lidar System 

Dr Ronald Fegley of the NOAA Wave Propagation Laboratory (Boulder, Colorado) is 
developing a dye lidar for installation at Point Barrow, Alaska (71° 20’ N, 156° 38’ W) as part 
of NOAA’s Geophysical Monitoring for Climatic Change program Starting in spring 1978, 
observations extending into the stratosphere are planned on a weekly basis, cloud cover permit- 
ting We intend to include such data as are made available in the SAM-II and SAGE ground 
truth data sets when this is appropriate However, the high probability of cloud cover at Point 
Barrow makes the probability of data capture within the time window of a SAM-II or SAGE 
overflight smaller than is acceptable for a primary ground truth sensor Moreover, the short 
wavelength (585 nm) of the lidar yields expected errors in measured particulate backscattering 
of greater than ±30% throughout the stratosphere under nonvolcanic background conditions 
(These errors result primarily from molecular density uncertainties, see also Section 2 2 1) If is 
for these reasons that an airborne lidar of longer wavelength (694 or 1060 nm) has been 
developed specifically for SAM-II and SAGE ground truth measurements 


2.3.4. Polar Nephelometer 

The Georgia Tech laser polar nephelometer (Grams, et al , 1975) measures the angular 
scattering pattern (often called the scattering phase function) of the light scattered out of a col- 
limated beam by aerosol particles This parameter is important for use in calculating the effect 
of airborne particles on the transfer of radiation in the atmosphere It can also be used to infer 
other physical properties of the aerosol particles, such as their size distribution or refractive 
index (e g , Grams, et al , 1974) 
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The nephelometer has been designed to operate in a pressurized aircraft cabin, using out- 
side air ducted through an airflow tube The sample volume is the intersection of a collimated 
source beam and the detector field of view within the airflow tube The source is a linearly 
polarized laser beam The optical system defines a collimated field of view (0 5° half-angle), a 
photomultiplier tube is located immediately behind an aperture in the focal plane of the objec- 
tive lens to measure the amount of light scattered from the laser beam A two-channel pulse 
counter, synchronized to the laser output, measures the photomultiplier pulse rate with the 
light beam both on and off The difference in these measured pulse rates is directly propor- 
tional to the intensity of the scattered light from the volume in which the laser beam and the 
detector field of view intersect 

The nephelometer is operated under the control of a microprocessor system The light- 
scattering measurements can be made at angles of 15°-165° from the direction of propagation of 
the laser beam Intermediate angles between these extremes are obtained by selecting the 
angular increments desired, pulses provided by digital circuits control a stepping motor that 
rotates the detector sequentially by preselected angular increments (usually 5°) The synchro- 
nous photon-counting system automatically begins to measure the scattered-light intensity 
immediately after rotation to the new angle has been completed 


2.3.5. In Situ Particle Sizing Devices 

2.3. 5.1. Quartz Crystal Microbalance Cascade Impactor 

The quartz crystal microbalance cascade impactor (QCM) is a multistage impactor that 
senses the mass of suspended particles as a function of particle size The particles are drawn 
into the sensor and separated aerodynamically into ten size intervals ranging from 0 05 /xm to 
25 /an in diameter (assuming spherical particles of mass density 2 gm citT 3 ) Figure 23 illus- 
trates the flow through the sensor The air velocity increases as it flows from one stage of the 
cascade to the next because of the decreasing diameters of the entrance jets Thus, the larger 
particles impact on the upper stages and the smaller ones are carried with the flow to the lower 
stages Table 8 lists the 50-percent cut points (diameter at which the particles have 50-percent 
impaction efficiency) for each of the ten stages Each impactor stage contains a piezoelectric 
crystal microbalance that senses the mass of the particles collected by a change in oscillator fre- 
quency between a reference crystal and the sensing crystal The time response is on the order 
of two seconds, so that temporal or spatial variations are resolved Laboratory analysis (scan- 
ning electron microscopy) can be performed on the samples for elemental composition and 
morphology 

This technique can provide data that are free of some of the ambiguities present in light- 
scattering and other remotely sensed data Information needed to expand the usefulness of 
satellite data for radiative transfer and associated modeling studies can be obtained 

NASA Langley owns two of these QCM instruments One of them, with an improved 
design for high-altitude sampling, was flown aboard the NCAR Sabreliner up to an altitude of 
13 km, it performed successfully in measuring the stratospheric aerosol size distribution over 
the Fairbanks, Alaska area in May 1977 and over the Sondrestrom, Greenland area in 
November 1978 as part of the SAM-II ground truth program It was also flown aboard the 
NCAR Queen Air over Guatemala, where, in February 1978, it successfully measured the size 
distribution of volcanic aerosols In addition, the instrumention was flown with the airborne 
lidar on the WFC P-3 during the Laramie practice comparative experiment in September 1978 
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SENSING CRYSTAL 


REFERENCE CRYSTAL 



FIGURE 23 ILLUSTRATION OF AIR FLOW THROUGH THE 10-STAGE QCM 
CASCADE IMPACTOR 


These experiments demonstrate the impactor’s suitability for airborne measurements and its 
ability to provide useful data on atmospheric aerosols Figure 24 shows the aircraft mounting 
for the microbalance Listed below are the specifications relevant to aircraft operation 


Cascade and Aerodynamic Housing 
External dimensions 

Weight 
Control Units 

External dimensions 

Weight 

Power Requirements 

Pump Motor 
Electronics 


5-inch-diameter cylinder 

22 in long with hemispherical end caps 

15 lbs 


19-inch rack mount 

7 in high 
12 in deep 

8 lbs 


28 Vdc, 2 amps 

115 Vac, 50-60 Hz, 1 amp 
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Table 8 


50%-EFFICIENCY POINTS FOR THE CELESCO MODEL C-1000 
QCM CASCADE IMP ACTOR 


Stage No 

50%-Cut-Point Diaijieter 
(Micrometers) 

1 

25 0 

2 

12 5 

3 

6 3 

4 

3 2 

5 

1 6 

6 

08 

7 

04 

8 

02 

9 

0 1 


0 05 


* 


Assuming spherical particles of mass density 2 g cm 


-3 


Major Equipment Subsystems 

Ten-stage cascade-sensing stack assembly with remote-controlled 
motor-driven inlet valve 

Sample-air pump 
Isokinetic air inlet 
Aerodynamic housing 

Control unit for 19-inch rack mounting (includes power supply, 
signal conditioning electronics, and printer) 

Interfacing cable 

The current plan is to transfer this high-altitude QCM to the NASA Ames U-2 and fly it 
as a SAGE correlative sensor at White Sands in February or March 1979 (See Experiment 2 in 
Table 2) 


2 . 3 . 5 . 2 , Other In Situ Particle-Sizing Devices 

It is possible that other particle-sizing devices may be employed on aircraft platforms dur- 
ing the SAGE ground truth program For example, a Climet single-particle optical counter was 
installed on the NCAR Sabreliner aircraft along with the Georgia Tech polar nephelometer and 
the NASA Langley quartz crystal microbalance during the November 1978 SAM II ground 
truth program in Sondrestrom, Greenland The principle of operation of this type of instru- 
ment 7 is similar to that of the University of Wyoming dustsonde (Section 2 3 1), i e , aerosol 


^Counters of this type are manufactured by such companies as Bausch and Lomb, 820 Linden Avenue Rochester New York 14625 
Climet Instruments Company 1240 Birchwood Drive Sunnyvale, California 94086 Particle Measuring Systems 1855 S 57th Court 
Boulder Colorado Roy to Instruments Inc 141 Jefferson Drive, Menlo Park California 94025 
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FIGURE 24 CONFIGURATION FOR EXTERNAL AIRCRAFT MOUNTING OF QUARTZ CRYSTAL 
MICROBALANCE 



particles flow through an illuminated volume, light scattered by an individual particle into a par- 
ticular solid angle is sensed photoelectrically, and the response pulses are classified according to 
their magnitude 

Other possible techniques that could be employed for in situ measurements include the 
use of devices for collecting particles by impaction and then analyzing the sample in the labora- 
tory to establish the particle size, shape, and composition by electron microscopy (e g , Ferry 
and Lem, 1974) In a similar fashion, Nuclepore filters with very small pore diameters might 
be used to collect particles to obtain samples for analysis by scanning electron microscopy (e g , 
Patterson and Gillette, 1977) Other techniques for sampling aerosol particles have been 
described by Grams and Rosen (1978) However, many of these techniques have not been 
employed in stratospheric sampling programs, and thus modifications of the sampling equip- 
ment would probably be required 


2.3.6. Noctilucent Cloud Sightings 

The International Noctilucent Cloud (NLC) Program, in cooperation with the World 
Meteorological Organization, records and collates observations of noctilucent clouds (particulate 
layers in the 73-95 km altitude region) Mr E J Truhlar of the Canadian Atmospheric 
Environment Service (AES) collects observations from 60 Canadian and 16 U S stations, and 
Dr D H McIntosh of the University of Edinburgh (UE) collects them for Western Europe 
We have made arrangements to receive tabulations of sightings (with about one month’s delay) 
from the AES and UE and will attempt to make similar arrangements with others Sample 
tabulations, relevant correspondence, and a description of the observation characteristics are 
given in Appendix C 

Resolution and Accuracy 

Bronshten and Grishiken (1975) estimate that the height accuracy of the "best" NLC 
measurements is better than 1 km (see also Appendix C) Presumably the "best" measure- 
ments are those made with aligned cameras or visually with theodolites Many sightings are 
made with the unaided eye, and their height accuracy is probably considerably worse than 1 km 
Uncertainties in the size distribution and shape of NLC particles, together with the uncalibrated 
brightness scale of observations, make conversion of sightings to 1 0-/um extinction coefficients 
extremely approximate at best The purpose of the NLC sightings will be to provide time and 
location data for comparison with any layers detected by SAGE in the 73-95 km region 
Because of the paucity of other particulate data m this height region, the NLC data are con- 
sidered an important source of ground truth information 


2.3.7. SAM II 

The Stratospheric Aerosol Measurement II (SAM II) was launched on the Nimbus 7 satel- 
lite in October 1978 SAM II is a scanning solar photometer that operates similarly to SAGE, 
however, it has only one optical channel--at wavelength 1 0 /am (the longest SAGE 
wavelength) SAM ll’s mission is to map vertical profiles of aerosol extinction by scanning 
through the atmosphere during spacecraft sunrise and sunset events Because of the Nimbus 7 
orbit these vertical profiles will be restricted to two latitude bands 64°-80° N and S (See, e g , 
Russell et al , 1978 ) 
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Since SAGE’s coverage is expected to extend from the equator to 75° N and S, there is 
some overlap in the coverage of the SAM II and SAGE sensors Moreover, the seasonal 
dependence of coverage indicates that there will be occurrences of near-coincidence in space 
and time between SAM II and SAGE scans When this occurs comparisons will be made— as a 
means of evaluating the mutual consistency of the two measurements In particular, near- 
coincidences are expected to occur in late spring and early fall of 1979 near Sondrestrom, 
Greenland SAGE Correlative Experiment 4 (see Table 2 and Section 3 2 4) aims to provide 
correlative ozone, aerosol, and other measurements at the time of one or more of these coin- 
cidences 


2.4. N0 2 and Multiconstituent Sensors 

2.4.1. Pepin Balloon-Borne Sunphotometer and Spectrometer 


The University of Wyoming has developed a balloon-borne sunphotometer and spectrom- 
eter system The photometer contains a seven-channel instrument, four channels of which are 
spectrally equivalent to the SAGE instrument The spectral wavelengths for the instrument are 

Channel # Wavelength (/tin) 

1 0 39 

2 0 43 

3 0 45 

4 0 60 

5 0 82 

6 0 93 

7 1 00 


The spectrometer scans the spectral interval between 0 37 and 1 1 /im with wavelength resolu- 
tion on the order of 0 05 /i.m During flight both instruments on the balloon gondola lock onto 
the sun and view it during the sunrise or sunset event 


Resolution and Accuracy 

The photometer and spectrometer have the capability to measure the intensity of the sun 
as viewed through the atmosphere with a precision of ±0 1% of that intensity Both instru- 
ments are calibrated during the extinction balloon flight by observing the sun at high elevation 
angles from the balloon platform and then extrapolating the observed signal to conditions out- 
side the atmosphere Observations are made to less than 0 01 airmass 

Both the photometer and spectrometer measuring systems make use of the total solar 
disk, rather than the partial disk observed by the SAGE scan system Since these systems are 
used at balloon altitude, not at spacecraft altitude, they achieve vertical resolution of the same 
order as SAGE, but do not require high pointing accuracy and are less sensitive to refraction 
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Conversion to Particulate and Ozone Extinction Profiles 

The measurements from the photometer and spectrometer systems will be converted to 
extinction profiles using the method outlined by Pepin (1970, 1977) The spectral intervals 
around the SAGE bands can be studied by using the spectrometer measurements made at spec- 
tral intervals near the SAGE observations 


2.4.2. Murcray Interferometer and Spectrometer 

NO; was first detected in the lower stratosphere by means of absorption features observed 
in infrared solar spectra, these features are associated with the v-$ band in the 1618 cm -1 region 
(Murcray et al , 1968, Goldman et al , 1970) The few altitude profiles currently available for 
NOj in this wavelength region have been derived from infrared solar spectra obtained by bal- 
loons, although some recent results have been reported for the same spectral region used by 
SAGE (Kerr et al , 1977, Goldman et al , 1978) 

In view of the small amount of data available concerning the altitude distribution of NO^, 
there is no accepted measuring technique that can be applied to obtain "ground truth" data for 
SAGE It is rather proposed that a balloon flight be performed to obtain a vertical profile of 
NO; using the infrared solar spectral technique for comparison with the profile obtained by 
SAGE The spectral data required for this comparison will be obtained by means of the 
University of Denver balloon-borne interferometer system, which consists of an interferometer 
capable of obtaining solar spectra with a resolution of 0 01 cm -1 covering the 1600 cm -1 region 
Previous experiments have shown that this resolution is more than adequate to determine an 
NO,> profile from sunset solar spectra The system includes a biaxial pointing control, a PCM 
telemetry system, an on-board digital magnetic tape-recording system, and power supplies for 
operating the various components The balloon flight will be performed as part of a series of 
flights for NASA and will obtain comparative data Many constituents in addition to N0 2 will 
also be measured 

NOi total column data have also been obtained from ground-based visible spectra with a 
technique initially proposed by Brewer (1973) and used extensively by Noxon By obtaining 
data on the ground from the scattered sunlight after sunset (solar zenith angles > 90° ) one 
can use an Umkehr technique to obtain some profile data The spectral resolution needed in 
this type of measurement is not high, the main requirement for the measurement is wide 
dynamic range 

A balloon-borne UV-visible grating spectrometer system suitable for collecting data of this 
sort will be incorporated into the P3 aircraft instrumentation The advantages of operating the 
system from the aircraft are twofold the increased latitude range over which comparison data 
can be obtained and the ability to obtain data when the weather precludes doing so from the 
ground The spectrometer system is a 0 5-meter Czerny-Turner system that is double-passed 
The radiation is interrupted by a tuning-fork chopper after the first pass and the radiation is 
synchronously detected, which results in a significant reduction in stray light The system is 
capable of a 0 3 A resolution, but for this application it will be degraded to -5 A Data will be 
recorded by means of a digital magnetic-tape recording system 
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2.4.3. LIMS Instrument Package 

The LIMS Instrument Package (LIP) consists of a number of instruments that will pro- 
vide measurements for validating the data obtained from the Limb Infrared Monitor of the 
Stratosphere (LIMS) experiments on Nimbus 7 This balloon-borne package is scheduled for 
two or three flights in support of the LIMS experiment The sensors on the LIP will measure 
ozone, water vapor, nitric acid, temperature, and pressure The instruments are (l)a Dasibi 
ozone sensor (Johnson Space Center), (2) an infrared emission radiometer (Belgium), (3) an 
absorbing filter sensor (NCAR), (4) a water vapor infrared radiometer (United Kingdom) and 
(5) an electrochemical concentration cell sonde and radiosonde (Wallops Flight Center) Max- 
imum altitude of the LIP is approximately 35 km All LIP data, as well as other ground truth 
data for LIMS, are expected to be stored at NOAA/NESS in Suitland, Maryland The sensor 
scientists are (1) Dr Don Robbins (JSC), (2) Dr Carlos Lippens (Belgium), (3) Dr Bruce 
Gandrud (NCAR), (4) Dr John Harries (UK) and (5) Mr Lawrence Rossi (WFC) Mr Ed 
Szajna and Mr Fred Witten (GSFC) are the LIP manager and engineer, respectively 

The possibility of making a LIP flight from Palestine, Texas in coincidence with a SAGE 
tangent scan is being explored with members of the European ad hoc ground truth group (see 
Table 3) In addition, a LIP flight scheduled for LIMS support in Cold Lake, Alberta, in Febru- 
ary 1979 may be sufficiently close in space and time to a SAGE tangent scan to provide useful 
correlative data 


2.4.4. Noxon Spectrometer 

As indicated above (Section 2 4 3) N0 2 was initially detected m the stratosphere by means 
of absorption features observed in infrared solar spectra obtained from a balloon N0 2 also has 
a strong absorption band in the visible region, Noxon has used this band to obtain total column 
data on stratospheric NO; with a ground-based spectrometer system The visible region has the 
advantage that spectral scans of the zenith sky made at solar zenith angles > 90° contain 
absorption features resulting from the path of solar radiation in the stratosphere Noxon’s sys- 
tem, designed to be portable, has been used by him to obtain data over a wide range of lati- 
tudes A description of his instrumentation, data reduction procedures, and the results he has 
achieved is contained in a recent series of articles (Noxon et al , 1978, Noxon, 1978) It is pro- 
posed that this system be utilized as part of the ground truth obervation program scheduled for 
late February 1979 at White Sands Proving Ground 


2.4.5. Schmeltekopf Chemiluminescent Sensor 

As indicated above, the amount of data currently available on the altitude distribution of 
N0 2 is limited All data obtained so far are based on spectroscopic techniques and have been 
secured by remote sensing Schmeltekopf has developed an instrument for obtaining in situ 
N0 2 data that makes use of a chemiluminescent NO sensor constructed for balloon-borne 
operation The unit first measures NO with the chemiluminescent technique and then, using a 
mercury lamp, converts the N0 2 present in the atmosphere to NO by photodissociation The 
NO concentration is then remeasured, the increase over the preceding measurement denotes 
the N0 2 concentration The system has not been flown yet, but will be flown as part of the 
SAGE ground truth plan if it is operational by then 
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2.4.6. LIMS 


The Limb Infrared Monitor of the Stratosphere (LIMS) experiment on the Nimbus 7 
satellite is being conducted to determine global-scale vertical distributions of temperature, as 
well as of several gases involved in the chemistry of the ozone in the stratosphere Profiles of 
ozone (O 3 ), nitrogen dioxide (NO;), nitric acid (HNO 3 ), water vapor (H 2 0), and temperature 
are determined with fine vertical resolution from the lower stratosphere (=10 km) to the lower 
mesophere (=65 km) These data are derived by inverting measured limb radiance profiles 
obtained by LIMS, an infrared multispectral scanning radiometer Measurements are made in 
each of six spectral regions one each in the 9 6 fim O 3 band, the 6 3 fun NOj band, the 6 2 
/am H 2 0 band, and the 11 3 /am HNO 3 band, and two in the 15 /am band of CO; 

This experiment is a follow-on to the successful Limb Radiance Inversion Radiometer 
(LRIR) experiment flown on Nimbus 6 to measure O3, H/), and temperature The LIMS 
instrument is identical to the LRIR in many respects, but with the essential difference that two 
detectors were added to the focal plane array and five parameters are being measured rather 
than three The horizon scan rate was also decreased from one degree per second to a quarter 
degree per second to provide improved signal-to-noise performance These changes facilitate 
the measurement of constituents with small signals (e g N 0 2 , HNO3) and allow extension of 
these measurements to lower and higher altitudes 

A programmed scanning mirror in the radiometer causes the field of view of the six detec- 
tors to make coincident vertical scans across the earth’s horizon The data from these scans are 
stored on tape for later transmission to the ground During data reduction the measured limb 
radiance profiles from the carbon dioxide channels are operated on by inversion algorithms to 
determine the vertical temperature distribution This inferred temperature profile, together 
with the radiance profiles in the other channels, is then used to infer the vertical distribution of 
the trace constituents 

Comparisons of the ozone and NO; profiles measured by SAGE with those measured by 
LIMS will be made for cases in which the two measurements are in sufficient proximity In 
addition, the LIMS Experiment Team plans to use SAGE and SAM-II data on aerosols and thin 
clouds to identify any possible effects on the LIMS data (Russell and Gille, 1978) 


2.4.7. SBUV/TOMS 

The Solar and Backscattered Ultraviolet (SBUV) instrument on the Nimbus 7 satellite 
should provide ozone profiles between 30 and 45 km with a vertical resolution of 3-5 km and a 
precision of 5-10% In this height region the time scale for ozone changes is expected to be 
several days or longer Therefore, the differing local measurement times of SAGE (sunrise and 
sunset) and SBUV (approximately local noon) are not expected to invalidate comparisons of 
observed ozone profiles It is planned that comparisons will be restricted to spatial separations 
of ±500 km Concurrent SAGE data on aerosols and on neutral density will be used to elim- 
inate some of the uncertainties inherent in the inversion procedures for the two experiments 

SAGE ozone profiles can also be compared against total ozone observations determined 
both by SBUV, which views an area approximately 200 km x 200 km, and by the Nimbus-7 
Total Ozone Monitoring System (TOMS), which views an area approximately 50 km x 50 km 
Total ozone observations will be compared during SAGE and SBUV/TOMS coincidences (as 
defined above) Moreover, the TOMS instrument will scan across the satellite track, thus 
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providing almost full global coverage for total ozone during a 24-hour period TOMS observa- 
tions may therefore also be used to indicate spatial variability at the time of the major ground 
truth experiments depicted in Table 2 


2.5. Temperature, Pressure, or Density Sensors 
2.5.1. Balloon-Borne Radiosonde 

The radiosonde 8 is a balloon-borne, battery-powered, meteorological instrument that 
automatically transmits to a ground receiving station radio signals relating to the pressure, tem- 
perature, and humidity of the air from the surface to a height of approximately 30 km Wind 
direction and speed data are calculated from the rise and horizontal drift of the radiosonde, 
which are measured by the ground receiving station (RAWIN Set AN/GMD) 

The radiosonde system is enclosed in a plastic and paper container and is approximately 
40 cm high, 17 cm wide, and 12 cm thick It weighs 392 grams without batteries Figure 25 
shows the flight configuration of this instrument when attached to a standard 1200-gram 
Weather Service balloon 

The precision of the radiosonde instrument is specified below 


Atmospheric Pressure 


Temperature 
Relative Humidity 


1060 mb to 50 mb 50 mb to 20 mb 20 mb to 2 mb 
±2 mb ±0 5 mb 0 25 mb 

±0 5°C 

Within ±5% at temperatures +40°C to 0°C 


Section 4 2 4 1 describes the radiosonde data, their management, distribution, and so forth 


2.5.2. Super Loki Datasonde 

The Super Loki datasonde is the standard rocket-launched meteorological payload used by 
the U S meteorological community It is used to measure vertical profiles of atmospheric tem- 
perature and wind between 70 and 20 km Figure 26 depicts this system in descent 
configuration suspended below its decelerator after launch by the Super Loki rocket This sys- 
tem is used in conjunction with a GMD set at 1680- MHz 

The precision of the Super Loki datasonde system is from 1 5 to 3 5 C for temperature 
data, ±3 msec -1 for the wind data, 1 5-3% for pressure and 3-5% for density 

Section 4 2 4 2 describes the data handling aspects of this instrument system 


^Vehicles and Sensors of ihe UPN 607 Applications Sounding Rocket Program October 1977 Preliminary 
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BALLOON 


TOTAL LENGTH 

~ 20 m 


PARACHUTE 


RADIOSONDE 


NOT TO SCALE 


FIGURE 25 RADIOSONDE FLIGHT CONFIGURATION 


5 1 






FIGURE 26 SUPER LOKI DATASONDE FLIGHT PROFILE AND GMD TRACKING SCENARIO 




2.5.3. NMC Global Data Net 


The National Meteorological Center (NMC) will supply SAGE with profiles of tempera- 
ture, pressure-surface heights, and density inferred from the NMC upper-air data network 
(principally balloonsondes, rocketsondes, and satellite sensors) Procedures for deriving these 
profiles and the data format are described in Section 4 2 4 4 
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3. MEASUREMENT SITES, SCHEDULES, AND LOGISTICS 


3.1. Aerosol Practice Comparative Experiment 

An aerosol Practice Comparative Experiment was held on September 27, 1978 in Laramie, 
Wyoming It included flights by the airborne lidar and a dustsonde, as well as ground-based 
lidar measurements at NCAR in Boulder, Colorado (see Section 2 3 3 2) Its purpose was to 
develop multisensor and multiplatform coordination procedures, and also to test data reduction 
and comparison techniques, as mentioned previously in Section 1 2 

Laramie, Wyoming, was chosen as the preferred site for the Practice Comparative Experi- 
ment because of its proximity to the University of Wyoming This was a cost-reducing factor 
that also simplified the logistics of the dustsonde balloon flight and maximized chances of 
instrument recovery Although the balloon-borne sunphotometer was unable to participate in 
the Practice Experiment, it is expected that both the sunphotometer and spectrometer will be 
available for the February 1979 ground truth experiments at White Sands 


3.2. Postlaunch Measurements 

3.2.1. Predicting SAGE Tangent Times and Locations 

The SAGE tangent times and locations are determined by the AEM-B orbit parameters 
(inclination, height, and shape) and the launch time and date (which determine the relative 
positions of SAGE, the sun and earth surface sites) Since none of these factors can be known 
precisely before launch, exact predictions of tangent times and locations can be made only after 
launch, when orbit parameters have been accurately determined Possible variations in orbit 
parameters are considerable because the Scout launch rocket is less controllable than some 
other vehicles, and AEM-B cannot make postinsertion orbit corrections Nevertheless, as a 
means of preparing for postlaunch contingencies, it is useful to predict tangent occurrences for 
nominal orbit parameters and launch times, as well as for probable variations in these factors 

Edwin F Harrison of NASA Langley Research Center has prepared such a set of 
parametric predictions Figure 27 shows examples of these results, as latitude-vs-longitude plots 
of all tangent locations for each of 12 months after launch, assuming a launch time of 1030 
EST, January 25, 1979, a circular orbit with a height of 600 km, and an inclination of 55° Also 
shown on each plot are circles of 500-km radius each, centered on Sondrestrom, Greenland, 
Wallops Island, Virginia, White Sands, New Mexico, and Natal, Brazil (The circles appear as 
ellipses because of the map projection ) These circles show at a glance the number of times that 
SAGE sunrise and sunset tangent scans will occur within 500 km of each site in a given month 

A more detailed look at this question of tangent opportunities is given by Table 9, which 
focuses on the first 60 days of the mission (a time of intense ground truth activity, see Table 2) 
and the sites of Wallops Island, White Sands, and Natal Table 9 considers a range of launch 
times, two launch dates, and two orbit shapes A list of prioritized requirements for SAGE 
tangent locations and times is given in Table 10 In view of Table 9, it appears that a launch 
time of about 1030 EST would yield the greatest probabilty of satisfying the requirements of 
Table 10 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS 

Calculations assume i = 55°, k = 600 km, launch time = 1 030 EST, 25 January 1 979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume i = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume i = 55°, k = 600 km, launch time ■ 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume i = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 



TANGENT LATITUDE, DEGREES 



FIGURE 27 


SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume i = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume t = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume i = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 


90 


ON 



• SUNRISE 
+ SUNSET 


90 1 20 150 180 210 240 

TANGENT LONGITUDE, DEGREES 
(h) AUGUST 1979 


360 


FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume i = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume i = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume i = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume i = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Continued) 

Calculations assume i = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses 
have 500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 
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FIGURE 27 SAGE SUNRISE AND SUNSET TANGENT LOCATIONS FOR SELECTED TIME PERIODS (Concluded) 

Calculations assume i = 55°, k = 600 km, launch time = 1030 EST, 25 January 1979 Ellipses have 
500-km radius and are centered at Sondrestrom, Wallops Island, White Sands, and Natal 



Table 9 


SAGE GROUND TRUTH OPPORTUNITIES AS A FUNCTION OF LAUNCH TIME 


Sage 

Launch 

(EST) 

OPPORTUNITIES WITHIN 500 km RADIUS 

Wallops Island 
(37 85 N, 75 48 W) 

White Sands 
(32 38 N, 106 48 W) 

Natal, Brazil 
(5 87 S, 35 32 W) 

Mission 

Days 

Distance, 

km 

Mission 

Days 

Distance, 

km 

Mission 

Days 

Distance, 

km 

0900+ 

40 3 

138 8 

13 9 

373 4 

102 * 

225 0 




174 * 

253 5 

47 2 * 

223 6 




41 4 * 

198 4 

55 7 

380 5 

093 0 f 

128 

464 0 

174 * 

443 9 

122* 

460 9 


19 3 * 

336 9 

43 4* 

447 0 

21 7 

489 5 

1000+ 

148 

388 7 

129 

439 2 

— 

— 


21 3 * 

482 0 

194 * 

4175 






45 4 * 

495 3 



1030 f 

16 8 

278 9 

179 

466 4 

— 

— 


22 2 * 

411 9 

21 4* 

369 1 



1 100+ 

188 

1586 

19 9 

379 6 

— 

— 


45 2 * 

432 1 






46 3 * 

460 0 





1130 t 

20 8 

36 1 

21 8 

340 5 

28 6 

419 3 


24 2 * 

435 5 

48 3 * 

370 1 




47 2 * 

383 0 






48 2 * 

471 1 





0900 f 

173 * 

376 3 

57 9 

269 9 

102* 

96 0 


183 * 

424 8 

58 9 

329 4 

187 

411 5 


37 3 * 

370 7 






38 3 * 

167 7 






59 8 

448 1 





093 0^ 

109 

432 3 

59 9 

364 1 

112* 

340 3 


11 8 

371 0 



197 

393 6 


183 * 

302 3 



20 7 

358 1 


193 * 

280 5 






39 3 * 

46 6 





looof 

128 

366 5 

194 * 

294 8 

21 7 

391 6 


13 8 

93 6 



22 7 

434 9 


148 

424 4 



48 2 * 

342 2 


20 3 * 

215 5 



57 2 * 

489 4 


21 3 * 

363 6 






41 3 * 

249 0 





1030f 

14 8 

321 6 

20 3 * 

190 3 

23 7 

483 6 


15 8 

158 8 

21 3 * 

407 5 

49 2 * 

364 2 


22 3 * 

377 1 

43 3 * 

456 0 

56 7 

193 7 


43 3 * 

456 9 






tnc = 55°, launch date = February 1, 1979 
*Sunset 

th = (600 X 750) km 
f h = (450 X 700) km 
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Table 9 (Concluded) 



OPPORTUNITIES WITHIN 500 km RADIUS 

Sage 

Wallops Island 

White Sands 

Natal, Brazil 

(37 85 N, 75 48 W) 

(32 38 N, 106 48 W) 

(5 87 S, 35 32 W) 

Launch 

(EST) 

Mission 

Distance, 

Mission 

Distance, 

Mission 

Distance, 

x Days 

km 

Days 

km 

Days 

km 

0900+ 

17 32 * 

331 1 

39 4 * 

413 2 



— 


60 8 

463 4 

40 4 * 

405 2 



0930 f 

128 

337 4 

149 

368 3 

11 2 * 

367 0 


13 9 

426 5 

174* 

104 7 

20 2 

305 3 






47 2 * 

341 1 

1000 f 

20 3 * 

88 1 

42 4* 

374 2 

— 

— 


41 3 * 

402 9 

43 4* 

202 1 



1030 f 

15 8 

427 3 

179 

200 8 

14 2 * 

466 1 


16 8 

143 0 

20 3 * 

167 0 

15 2 * 

424 7 






24 2 * 

256 3 






58 7 

324 2 

1 100+ 

23 2 * 

175 3 

45 3 * 

462 8 

— 

— 


44 3 * 

240 0 

46 3 * 

127 8 




45 3 * 

395 3 





0900 f 

99 

370 8 

174* 

3188 

— 

— 

0930^ 



39 4* 

458 3 



128 

499 4 

174* 

388 4 

11 2 * 

426 6 




184* 

483 3 

20 7 

2164 






47 2 * 

303 8 

-++— 

O 

o 

o 





55 7 

493 9 

128 

374 1 

149 

423 6 

— 

— 


20 3 * 

367 3 

42 4 * 

343 7 






43 4* 

3102 



1030* 

15 8 

437 0 

179 

420 2 

15 2 * 

375 0 


168 

450 7 

20 3 * 

434 0 

23 7 

446 3 




21 4* 

394 5 

24 7 

363 3 






50 2 * 

350 3 

1100* 





58 7 

272 2 

15 8 

450 0 

179 

4147 

— 

— 


16 8 

404 4 

45 3 * 

365 5 




23 2 * 

188 8 

46 3 * 

164 6 




24 2 * 

429 8 






44 3 * 

313 8 






100 km= 

1 1 l°lat 

1 00 km= 

1 07°lat 

100 km= 

0 91°lat 


= 

0 87°long 


0 86°long 

= 

0 86°long 


inc = 55°, h = 600 km 
♦Sunset 

| Launch date = February 1, 1979 
'Launch date = January 25, 1979 
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Table 10 


PRIORITY OF SAGE SATELLITE COINCIDENCE SITES 
FOR GROUND TRUTH MEASUREMENTS 


1st 


2nd 


3rd 

4th 


No unocculted period during first two months to avoid interference with initial ground 
truth experiments (listed below) 

Sunrise and sunset satellite profiles occurring within a 5-day period at White Sands, NM 
(32°23’ N, 106°29’ W) during mid-February to early March 1979 

— or — 


Sunrise profile at White Sands 


— or — 


Sunset profile at White Sands 

Sunrise and/or sunset satellite profile at Wallops Island, Virginia (37°51’ N, 75°29’ W) 
during mid-February to early March 1979 

Sunrise or sunset satellite profile at Natal, Brazil (5°52’ S, 35°19’ W) during mid-to-late 
March 1979 


The results shown in Figure 27 and in Tables 9 and 10 have been taken into account in 
choosing the experiment sites and schedule shown in Table 2 Definitive tangent locations and 
times will be recomputed after launch, when more precise orbit data are available These results 
will be made available to all correlative sensor scientists, including those on the ad hoc ground 
truth teams and others outside the United States 


3.2.2. Required Proximity of SAGE and Correlative Measurements 

The required proximity of SAGE tangent scans and correlative measurements is deter- 
mined by the typical variability of the measured constituents On the scales of interest (several 
hundred km and several hours) this variability is difficult to assess because of the very small 
number of appropriate measurements that have been made During nonvolcanic conditions lidar 
measurements of stratospheric aerosols have frequently shown that, within the layer of max- 
imum aerosol content (-16-24 km), profile shape and magnitude were typically preserved 
throughout a night of observations, when stratospheric wind velocities were about 10 msec -1 
Assuming an observation time of about eight hours, this converts to a spatially uniform region 
of about 300 km or more Furthermore, occasional comparisons have been drawn between stra- 
tospheric lidar measurements made on the same night at locations separated by 1500-2400 km 
(17-27° of longitude and 2-5° of latitude, the locations included Menlo Park, California, 
Laramie, Wyoming, Boulder, Colorado, and Kansas City, Missouri ) On these occasions, which 
were 8 or more months after any noticeable volcanic injections, approximate spatial uniformity 
was also observed 

In view of these results, it appears that a proximity of ±500 km and ±2 hr between 
SAGE and correlative aerosol measurements would be satisfactory during times not appreciably 
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perturbed by volcanic activity Conditions can be expected to be more variable during the very 
interesting period just after a volcanic injection However, flights by the P-3 lidar will be made 
to document any spatial variations in aerosol structure between the SAGE tangent location and 
the site of any balloon flights 

Diurnal changes in ozone concentration caused by changes in solar irradiance should be 
primarily confined to altitudes above 50 km, which are above the range of SAGE measure- 
ments Between 50 and 25 km, ozone changes are related both to temperature- induced chemi- 
cal changes and to transport variations Below 25 km, where ozone variations are related 
almost entirely to transport variations, extremely large ozone variations occur, the vertical dis- 
tribution found on one day may bear little similarity to that of the previous day These low- 
altitude changes are primarily responsible for the variations in total ozone from day to day and, 
although these variations are difficult to predict, a high correlation is found between total ozone 
and changes of the 100-mb height field The day-to-day change in maximum ozone concentra- 
tion is about 20%, and the altitude of the peak may change by as much as 4 km 

In view of this temporal variability, correlative ozone measurements should be taken as 
close in time to the SAGE measurements as possible, probably within ±2 hours of the tangent 
scan It is difficult to translate the temporal variations described above into their spatial counter- 
parts However, the spatial variations can be expected to be meteorology-dependent and possi- 
bly significant For this reason the LIMS Correlative Measurements Plan specifies a desirable 
proximity between LIMS scans and correlative sensors of ±0 5° latitude and ±1° longitude, 
with maximum separations of ±2° latitude and longitude At the SAGE ground truth sites, 1° 
of latitude or longitude is roughly 100 km (to within ±14%, except at Sondrestrom, where 1° 
of longitude is only about 40 km) Thus, the LIMS plan specifies maximum separations of 
about ±200 km From Figure 27 and Table 9 it can be seen that, within any given month, the 
probability of a SAGE tangent scan’s falling within 200 km of a ground truth site is very small 
Thus, it would be unrealistic to specify such a small coincidence requirement for SAGE correla- 
tive measurements Nevertheless, because of the variabilities mentioned above, it is important 
that each SAGE correlative experiment attempt to characterize the local spatial variability near 
the time of the experiment This will be done by making ozone profile measurements several 
days and hours before and after SAGE tangent scans, as well as during the scans (See Sections 
3 2 4 and 3 2 5) 


3.2.3. Cluster Concept 

Because of the need to measure multiple constituents (hence, to use multiple sensors) 
and the desire to take advantage of existing capabilities, the SAGE Ground Truth Plan was 
built around measurement facilities that formed natural clusters These clusters are described in 
more detail in Table 3 (Section 1 2) Within some of the clusters simultaneous aerosol or ozone 
measurements can be made from different sites Such simultaneous measurements will be use- 
ful in characterizing spatial variability in the region of SAGE tangent scans 
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3.2.4. Northern Hemisphere Measurements 


3.2.4. 1. Experiment 1— Wallops Island/Hampton, April 1979 

The sensors to be included in this experiment have been listed in Table 2 Figure 28 
shows a typical measurement sequence for the correlative sensors that are launched within a 
few hours of the SAGE tangent scan 


3. 2. 4.2. Experiment 2-Holioman/White Sands, March 1979 

Sensors are listed in Table 2 and more information is presented in Table 11 A typical 
time sequence both for sunrise and for sunset scans is shown in Figure 29 


3. 2. 4. 3. Experiment 4— Sondrestrom or Poker Fiats, Spring or Summer 1979 

Sensors are listed in Table 2 A typical measurement sequence both for SAGE sunrise and 
for SAGE sunset scans is shown in Figure 30 Note that, because of its retrograde orbit, SAM 
II views a sunset while SAGE and earthbound observers view a sunrise, conversely, SAM II 
views a sunrise while SAGE and earthbound observers view a sunset 


3. 2. 4. 4. Experiment 5— Wallops Island/Hampton, Fall 1979 

Sensors are listed in Table 2 The measurement sequence is as shown in Figure 28 


3. 2. 4.5. Experiment 6— Wallops Island/Hampton, Winter 1979-80 

Sensors are listed in Table 2 The measurement sequence is as shown in Figure 28 


3. 2. 4.6. Other Northern Hemisphere Measurements 

It is expected that other northern hemisphere measurements will be scheduled by the 
European and Japanese ad hoc ground truth groups (See Section 2 1, Tables 4 and 5, and Fig- 
ure 3 ) The SAGE Experiment Team will furnish predictions of SAGE tangent locations and 
times to these groups The correlative sensor scientists will then schedule observations on the 
basis of these opportunities and their own constraints 


3.2.5. Southern Hemisphere Measurements 

3. 2.5.1. Experiment 3— Foraleza/Natal, April 1979 

Sensors are listed in Table 2 The measurement sequence will be as shown in Figure 28, 
with the addition of a dustsonde sequence as in Figure 29 
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FIGURE 28 TYPICAL CORRELATIVE MEASUREMENT SEQUENCE FOR SAGE SUNRISE OR 
SUNSET PROFILE AT WALLOPS ISLAND 



Table 1 1 


PRELIMINARY DESCRIPTION OF SAGE GROUND TRUTH 
EXPERIMENTS AT HOLLOMAN AFB, NM 


Measurement 

Measurement 

Coordinator(s) 

Flight 

Schedule 

Flight 

Duration 

Maximum 

Altitude 

Maximum 

Diameter 

Load Train 
Length 

Transmitter 

and 

Payload 

Weight 

Airborne 

Lidar 

W H Fuller 
(NASA-LaRC) 

5 fits 

(1 flt/day) 

4 hrs 

16,000 ft 




Dustsonde 

D T Hoffman 
(U of Wyoming) 
J M Rosen 
(U of Wyoming) 

1 fit 

3 hrs 

100,000 ft 

50 ft 

320 ft 

20 lbs 

Balloon 

Sunphotometer 

T J Pepin 
(U of Wyoming) 

1 fit 

9 hrs 

100,000 ft 

60 ft 

20 ft 

150 lbs 

Rocket 

Ozonesonde 

T Perry 
D Bruton 
(NASA-WFC) 

2 fits 

2 hrs 

250,000 ft 

16 ft 

18 ft 

2 lbs 

Balloon 

Ozonesonde 

T Perry 
D Bruton 
(NASA-WFC) 

4 fits 

2 hrs 

100,000 ft 

9 ft 

50 ft 

4 lbs 

Balloon 

Interferometer 

D G Murcray 
(U of Denver) 

1 fit 

5 hrs 

120,000 ft 

275 ft 

200 ft 

1400 lbs 

U-2 

SO sub 3 $, NO, T 
Sensor 

M Lowenstein 
T Starr 
(NASA-ARC) 







U-2 Wire 
Impactor 

N Farlow, 

G Ferry 
(NASA-ARC) 

2 fits 

6 hrs 

70,000 ft 




U-2 Quartz 

Crystal 

Microbalance 

M P McCormick 
(NASA-LaRC) 

D C Woods 
(NASA-LaRC) 

2 fits 

6 hrs 

70,000 ft 




U-2 Multifilter 
Package 

A Lazrus 
(NCAR) 

2 fits 

6 hrs 

70,000 ft 




U-2 Glass 
Collector 

R Charlson 
(U of Washington) 

2 fits 

6 hrs 

70,000 ft 
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FIGURE 29 TYPICAL CORRELATIVE MEASUREMENT SEQUENCE FOR SAGE SUNRISE AND 
SUNSET PROFILE AT HOLLOMAN-WHITE SANDS IN FEBRUARY 1979 
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FIGURE 30 TYPICAL CORRELATIVE MEASUREMENT SEQUENCE FOR SAGE SUNRISE AND 
SAM-II SUNSET PROFILE AND FOR SAGE SUNSET PROFILE AT SONDERSTROM 
IN SUMMER 1979 


75 



3.2. 5. 2. Other Southern Hemisphere Measurements 

It is hoped that other southern hemisphere measurements will be scheduled by CSIRO in 
Australia and INPE in Brazil (See Appendix B, Table B-l ) In particular, Dr John Gras of 
CSIRO is planning balloon-borne measurements of aerosol concentration for several size 
cutoffs 
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4. DATA PROCESSING 


Procedures for reducing ground truth data to formats that are readily comparable to SAGE 
extinction and number profiles will be established and tested prior to launch, so that data com- 
parisons can be carried out in the shortest time possible 


4.1. Ancillary Data 

Both the lidar and sunphotometer data require a molecular density profile for separation 
of gaseous from particulate optical coefficients In addition, the lidar data (especially ruby data) 
require an estimated ozone and particulate extinction profile A two-step process will be used to 
incorporate these ancillary data into the data analysis 


4.1.1. Model Atmospheres 

Prior to launch, model profiles of atmospheric density and of ozone and particulate extinc- 
tion (at the lidar wavelengths) will be developed for the location and month of each ground 
truth site (see Table 2) These profiles will be stored on cards or another medium that can be 
computer-read as input to the lidar and sunphotometer data reduction algorithms (see Section 
4 2) They will be used for initial reduction of the lidar and balloon photometer data 


4.1.2. Measured Atmospheres 

Each ground truth experiment will include at least one measurement of the local tempera- 
ture and pressure profiles (These will be provided by standard radiosonde, by radiosonde pack- 
ages on the ozone balloonsondes, by the rocket-launched datasondes, and by the dustsondes, 
which measure temperature and pressure in addition to the aerosol data See also Section 
4 2 4 ) These profiles will be converted to density profiles and stored in the same format as the 
model density profiles, so that they can be readily substituted for the model profiles in the data 
reduction In addition, as a routine part of SAGE data reduction, density profiles for the time 
and location of each SAGE scan will be derivYd from the SAGE data These SAGE density 
profiles will be extracted from the SAGE data tapes and stored in the same format as the model 
and sonde-measured density profiles Model, sonde-measured, and SAGE-measured density 
profiles will be plotted for each ground truth experiment when available In this manner 
differences in density profiles can be highlighted and, by using each profile sequentially in the 
data analysis, effects of density differences on derived constituent profiles (both SAGE and 
correlative) can be explored and understood 

Carefully chosen model ozone profiles are probably adequate for lidar data reduction (cer- 
tainly for Nd lidar), nevertheless, measured ozone profiles near the time and location of lidar 
flights will be compared to the model profiles and, if necessary, substituted Likewise, particu- 
late extinction profiles (at the lidar wavelength) derived from the lidar and SAGE measure- 
ments will be compared to the model profiles and iteratively substituted whenever appropriate 
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4.2. Production of Data Products for Comparison to SAGE Data 


4.2.1 Ozone Data 

4.2. 1.1. Dobson Spectrophotometer 

Dobson spectrophotometer observations will normally be made three times daily (morn- 
ing, near local noon, and afternoon) The observational procedures that will be followed are 
outlined in the Weather Service Observer’s Manual 9 The Dobson dial readings and other per- 
tinent data are entered by hand on the NOAA observation form, NWS form B-35b (3-73) 
shown in Figure 31a Instructions regarding the use of this form are provided on the reverse 
side (Figure 31b) The total ozone overburden above each site is derived by applying the 
appropriate correction factors and the air mass data pertinent to that site Sample results 
obtained by both manual and computer data processing are shown in Figures 32 and 33 


4.2. 1.2. Canterbury Spectrometer 

The Canterbury spectrometer will conduct basically the same type of observations, with 
similar data output and data delivery schedules as for the Dobson spectrophotometer 


4. 2. 1.3 Balloon ECC Ozonesonde 

The data measured as the ECC ozonesonde rises consist of a vertical profile of ozone 
However, since the ozonesonde is attached to a standard radiosonde, atmospheric pressure, 
temperature, and relative humidity data are also measured (See Section 4 2 4) The ozone data 
are fed to the telemetry unit in the radiosonde and transmitted in flight to an AN/GMD, where 
they are automatically recorded along with the pressure, temperature, and humidity data on a 
TMQ-5 recorder The GMD antenna angles are recorded on a control recorder Figure 34 dep- 
icts the acquisition and flow of ozone data These are subsequently digitized and inputted to the 
Wallops ozone computer program, where they are merged with associated calibration data The 
calibration data are obtained from the ozonesonde’s preflight calibration process and the 
radiosonde’s preflight baseline checks The ozone data are then outputted as tables (Figure 35) 
and charts (Figure 36) Typically, the final data for flights launched from Wallops will be avail- 
able within 1-2 weeks Final data from launches away from Wallops (allowing for transmittal 
from the field site to Wallops) will be available within 3-4 weeks These data will also be pro- 
vided on magnetic tape Standard card and plotting formats for comparisons of SAGE and 
correlative profile data are described in Section 4 2 5 


4.2. 1.4. Super Loki Optical and Super Areas Chemiluminescent 
Ozone Payloads 

Both these systems are launched on rockets and are ejected at apogee to descend on a 
decelerator The ozone measurements are made during the descent phase Figure 37 shows the 
data flow for these systems The measured data are transmitted to an AN/GMD which is 
modified to allow it to receive and record an 8-bit PCM telemetry signal from these payloads 
The frequency is the standard AN/GMD 1680 MHz In addition, to obtain accurate position 
data of the ozone sensor during descent, radar tracking is required All the rocket systems 


^Observer’s Manual, Dobson Ozone Spectrometer revised November I 1972 U S Department of Commerce National Weather Ser- 
vice 
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The intapMM of VS Form B — 33b is dictated by rnmporrr pio gyamm mg aa4 poach card ftfMWWM Observational 
data aha aid b« entered only to the revs marked with aa asterisk, nod should b« racocdad ta successive columns, accept 
that blank columns lbs aid be used to aaparaca data obuuoad oo differ esc day a 

2. BNTtYOPDATA 


a. Va t l aa • Enter (he name of the nation 
k. tesfresmet Naa i har - Record (ha serial c 


r of the Dobaon Omae S p e cti opbmom eier employed. 


c. Observations - Enter (he s pectrophotometer dial feadmgs to (he acareat tenth ot a decree Observational procedures are 
eetlined ia the Observers' Handbook Dobson Ozone Spectrophotometer 

4 XododDete 


(1) Type ef Observation • Indicate the type of obsm 


Wove lengths Used 


ADDSGQP 

CDDSGQP 


■ mode, acmtd mg to the code m (he following coble 
Light t eas e s 
Direct Sue, Using GQP 
Danes Sm, Using GQP 


ADZC A and D dandy Zenirh 

CCZC C and C' Cloudy Zenith 
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Code Wavelengths Used Light Sonrce 
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Motes 
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C - Cleat H - Hazy VH - Very Hazy 

(3) Description of the Cleody Zenith - Then an observation is made on the cloudy zenith, indicate the cloud height, 
thickness and texture m the vicinity of the zenith according to the code in the following table 

Clewd Height Owed Thickness Clasd Texture 


Space ia provided under "Notes ' for recording additional information pertaining to the 
indicate the presence of fog, amok# etc 


P - Patchy 

accuracy of an observation, e g , 


L Obi error's Initials - The observer sirould until the obaermcsooai data (or which he is responsible 
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• TOTAL ozone REDUCTION PROGRAM 
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FIGURE 35 EXAMPLE OF COMPUTERIZED BALLOON OZONE DATA PRINTOUT 
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FIGURE 37 


DATA FLOW FOR SUPER LOKI OPTICAL OR SUPER ARCAS CHEMILUMINESCENT 
OZONE PAYLOAD 








planned for the SAGE ground truth program will be launched from sites at which support is 
provided from both AN/GMD and radar systems The ozone data are recorded on an analog 
tape Wallops will provide an analog magnetic recording system for recording the ozone data to 
interface with the AN/GMD at other sites The radar installations will record radar data on a 
digital magnetic tape Both the ozone and radar magnetic tapes, as well as the associated paper 
records, will be forwarded to Wallops along with the data from the supporting balloon 
radiosondes (Section 4 2 4), the ECC ozonesondes, and the Super Loki datasondes The rocket 
ozonesonde data are digitized at Wallops The rocket ozonesonde and radar data, along with the 
lower-level balloon radiosonde, ECC ozonesonde, and the upper-level Super Loki Datasonde 
data (temperature/winds), are then forwarded to Goddard Space Flight Center These data and 
the calibration data associated with the respective rocket ozonesonde are then merged and 
reduced to form the final data product In the case of the Super Loki Optical Ozonesonde, the 
data are reduced by computer and plotter The final products are an ozone data table, as shown 
in Figure 38, and plots of ozone density and the ozone mixing ratio, as shown in Figures 38 
and 39, respectively In the case of the Super Areas chemiluminescent ozonesonde, the data 
reduction is performed with a desk-top minicomputer and the final product is an ozone data 
table as shown in Figure 40 

After reduction, the data for both types of rocket system are forwarded to the Wallops 
ASRP office for packaging and distribution Typically, the data will be ready for distribution 
within eight weeks following the sounding 

Standard card and plotting formats for comparisons of SAGE and correlative profile data 
are described in Section 4 2 5 


4.2.2. Aerosol Data 
4.2.2. 1. Airborne Lidar 

The lidar data will be reduced to vertical profiles of the particulate backscattenng 
coefficient (at the lidar wavelength) by using data reduction techniques similar to those rou- 
tinely employed by SRI and NASA Langley in many previous measurements (See, e g Russell 
et al , 1976a, b) The reduction algorithm will automatically compute error bars that include 
uncertainties in (1) signal measurement, (2) density estimation, (3) transmission estimation, 
and (4) normalization (See Figure 41 ) In addition, scattering ratio profiles will be normalized 
to make the minimum scattering ratio equal the value expected on the basis of previous dust- 
sonde measurements and optical models (e g , Russell et al , 1976b), rather than the value of 
unity that has customarily been assumed in the past This procedure has the effect of sym- 
metrizing the expected normalization error and reducing it by about half 

Lidar-measured particulate backscattenng coefficients can be converted to particulate 1 0 
/i.m extinction coefficients by using an assumed refractive index and particle size distribution 
Figure 42 shows the dependence of the conversion ratio on optical model properties (cf Figure 
13 and accompanying discussion) (The size distributions and compositions shown have been 
derived from measurements by various investigators --eg Hofmann et al , 1975, Pinnick et 
al , 1976, Shettle and Fenn, 1976, Deirmendjian, 1969, Toon and Pollack, 1976, Harris and 
Rosen, 1976, Swissler and Harris, 1976) In a given lidar measurement the optical model can in 
general only be estimated on the basis of previous or simultaneous (e g dustsonde) measure- 
ments Numerous dustsonde measurements have shown a preferred height dependence for the 
channel ratio, N 1S /N (defined in Section 2 3 1) Specifically, for nonvolcanic conditions, 
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ROCKET OZONE DATA 


_Location Wallops Island 
p.n«nm»nnt ^- Krueger 

—Payload No 209 


_Rocket Total Ozone Above _ 1 

Equals • 

.Balloon Residual Ozone Below 13 » 5 | 

Equals 103 

_Total Ozone = • 31 7 

.293 3 1S40 

.Dobson Total Ozone- .296 3 1739 


0 3 

Mixing Partial 

Ax/ Ah Probable x(h) E(h) Ratio Pressure 

Ht km atm-cm/km Error % atm-cm motfcm 1 l/gm/gm ymb 




Air Temp, Pressure Oensity, Oala Source 1715Z at Wallops Island, VA, for 19-23 k m. , used 


the balloonsonde for 11/18/76 at 17262 at '.lallons 
Island, VA 


FIGURE 38 EXAMPLE OF ROCKET OZONE DATA 
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(a) DENSITY (MOLECULES/cm 3 ) 



FIGURE 39 EXAMPLES OF OZONE CONCENTRATION PROFILE MEASURED BY OPTICAL 
ROCKETSONDE 
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Poker Flat, Alaska September 25, 1976 0756Z 
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FIGURE 40 EXAMPLE OF SUPER ARCAS CHEMILUMINESCENT OZONE DATA 
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FIGURE 41 AIRBORNE LIDAR DATA-PROCESSING FLOW 
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E(1 0 pml/BlX) 



FIGURE 42 EXTINCTION-TO-BACKSCATTER RATIOS, E(1 Opiml/BU), COMPUTED FOR RUBY 
AND Nd WAVELENGTHS X, SILICATE AND AQUEOUS SULFURIC ACID 
COMPOSITIONS, AND A WIDE RANGE OF SIZE DISTRIBUTIONS 

Size distributions are grouped according to their value of the integral number ratio 
N is/N 2 s See Table 6 and text for discussion of observed regions and frequency 
of N js/N 25 values 
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FIGURE 43 DATA FLOW FOR BALLOON-BORNE RADIOSONDE PAYLOAD 




N 15/N 25 values within the major aerosol mixing-ratio peak are usually between 3 and 5 (Hof- 
mann et al , 1975) Thus the uncertainty in conversion ratio is given by the vertical spread of 
data points in Figure 42 above the appropriate range of N )5 /N 25 values The mean and stan- 
dard deviation of appropriate subsets of conversion ratios is shown in Table 12 

Figure 42 and Table 12 show that when recent or nearby dustsonde measurements, or 
other measurements, indicate that N ]S /N 25 falls in the range 3-5, the uncertainty (standard 
deviation) in the conversion ratio is about ±10% for both ruby (X=0 694 /im) and 
Nd(X=l 06 /iin) measurements The slightly larger uncertainty for ruby measurements shown 
in Table 12 arises from uncertainty in particle composition (aqueous sulfuric acid or silicate), to 
which the Nd conversion ratio is not so sensitive (presumably because the Nd wavelength is 
close to 1 0 yum) If the particle composition can be ascertained by some other 
measurement(s), then the uncertainty in the ruby conversion ratio is reduced slightly — to about 
±8%, which is about equal to the uncertainty in the Nd ratio 

The derived profiles of particulate backscattering coefficient will be converted to profiles 
of 1 0-yxm extinction coefficient by using an appropriate conversion factor from Figure 42 
Extinction error bars will also be computed, based on the particulate backscattering error bars 
and the uncertainty in the conversion factor 

Both the backscattering and extinction profiles will be plotted on standard scales and 
stored on standard-format punched cards, as described in Section 4 2 5 

The lidar-measured backscattering coefficient profile will be superimposed on the SAGE 
extinction coefficient profile to derive a cross-wavelength extinction-to-backscatter ratio (which 
may be height-dependent) Conversion ratios derived in this manner will be compared with the 
values in Figure 42 to facilitate the selection of appropriate optical models and subsequent 
aspects of data validation and reduction 


4. 2. 2. 2. Dustsondes 

The dustsonde data will be reduced to vertical profiles of N | 5 and N 15/N 25 by means of 
the data reduction techniques the University of Wyoming has routinely employed for many 
years (N x is the number of particles with radius > xyum, cf Figure 12) The derived profiles 
of N 15 and N 15/N 25 will be converted to profiles of 1 0-jim extinction coefficient by using an 
appropriate conversion factor from Figure 13 or 14 Extinction error bars, based on the uncer- 
tainties in N 15/N 25 and in the conversion ratio, will also be computed 

Number and extinction profiles derived from the balloon data will be plotted on the stan- 
dard scales described in Section 4 2 5 The extinction will be compared directly with the 
corresponding Ildar and SAGE results In addition, the particle number profile will be used to 
derive an extinction-to-number ratio, possibly height-dependent, for comparison with Figure 
13-and will be directly compared with the particle number profile derived from the SAGE 
extinction profile 


4.2. 2. 3. Balloon-borne Sunphotometer 

The sunphotometer data, using algorithms now being developed by T Pepin at the 
University of Wyoming, will be reduced to vertical profiles of particulate 1 0 j/.m extinction 
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Table 12 


MEAN AND STANDARD DEVIATION OF EXT1NCTION-TO-BACKSCATTER RATIO, 
E(1 Ojum)/B(X), FOR TWO LIDAR WAVELENGTHS X AND 
VARIOUS GROUPS OF SIZE DISTRIBUTIONS 



Lidar 

Ruby (X=0 6943 urn) 



Nd-YAG (X= 

06 jum) 


Composition 

N 15 /N 25 

1-2 

2-3 

3-5 

3 5-4 5 

5-10 

10-20 

>20 

1-2 

2-3 

3-5 

3 5-4 5 

5-10 

10-20 

>20 

75%H 2 S0 4 -25%H 2 0 

Mean (sr) 
o/Mean 

30 

33% 

38 

6% 

35 

8% 

35 

5% 

23 

21% 

14 

24% 

8 

25% 

44 

46% 

70 

9% 

63 

8% 

63 

3% 

43 

18% 

28 

14% 

21 

13% 

Silicate 

Mean (sr) 

23 

29 

31 

32 

23 

15 

9 

32 

61 

61 

61 

45 

31 

23 

a/Mean 

22% 

13% 

8% 

6% 

20% 

24% 

12% 

51% 

22% 

10% 

6% 

18% 

12% 

5% 

Both 

Mean (sr) 

26 

34 

33 

33 

23 

14 

8 

38 

65 

62 

62 

44 

29 

22 

a/Mean 

31% 

16% 

10% 

7% 

21% 

24% 

20% 

50% 

17% 

9% 

5% 

18% 

14% 

12% 




These results will be plotted on the standard scales and punched in the standard format Error 
bars will be provided at 1-km intervals 


4. 2. 2. 4. Polar Nephelometer 

The nephelometer records data on the scattering phase function for an aerosol Measure- 
ments of the light scattered from the laser beam at the wavelength (633 nm) of the helium- 
neon laser used in the instrument are recorded in 5-degree steps at angles between 15 and 165 
degrees When the nephelometer is used on an aircraft platform such as the NCAR Sabreliner, 
phase function measurements can be obtained at some three to five altitudes along SAGE 
tangent paths For data validation, analysis can be made of other in situ data obtained simul- 
taneously on the same aircraft with the quartz crystal microbalance, single-particle optical 
counter, and particle collection devices to determine particle size distribution, shape, and com- 
position In this way an effective value of the complex refractive index of the particles can be 
established by means of least-squares curve-fitting techniques similar to those described by 
Grams et at (1974) These results, together with the aerosol number density profiles obtained 
by the dustsonde, the backscatternng profiles obtained with the airborne lidar, and the extinc- 
tion profiles obtained with the SAM II and SAGE sensors, will be used to improve the optical 
model 

In addition, the intercomparisons among all the above measurements will allow inferences 
to be drawn regarding the radiative properties of the stratospheric aerosol layer These results, 
combined with the SAM II and SAGE data on the spatial and seasonal variability of the stratos- 
pheric aerosol layer, can then be used by the experiment team members in a variety of SAM II 
and SAGE data-use investigations Of particular interest for these investigations will be the 
so-called asymmetry factor, single-scattering albedo, and extinction coefficient for use in climate 
theories involving two-stream radiative transfer approximations (see, e g , Chylek and Coakley, 
1974) Estimates of these parameters will be based on analysis of the scattering data (obtained 
with the polar nephelometer and the lidar) and the extinction coefficients measured by the 
SAGE instrument 


4.2.3. NO 2 and Multiconstituent Data 

4.2.4. Temperature, Pressure, and Density Data 
4.2.4. 1. Balloon-borne Radiosonde 

The data measured as the radiosonde balloon rises consist of vertical profiles of atmos- 
pheric pressure, temperature, and relative humidity The data are transmitted in flight to a stan- 
dard ground-based AN/GMD (see Figure 26), 10 where they are automatically recorded on a 
TMQ-5 recorder The GMD antenna azimuth and elevation angles are recorded on a control 
recorder At sites without computer support an observer transcribes the TMQ-5 and control 
recorder’s chart data onto standard plot forms Measurements of pressure are made in millibars, 
temperature in Celsius, and moisture in percent of relative humidity Wind direction and speed 
are also determined in this process Figure 43 illustrates the flow of the radiosonde data Fig- 
ures 44, 45, and 46 comprise samples of the standard NOAA Weather Service Charts (Forms 
MF3-31A, B, C) used in compiling the basic raw data, while Figure 47 is a computerized reduc- 
tion of the data These data, whether from the manual or the computerized reduction, are 

'^Federal Meteorological Handbook No 3 Radiosonde Observations US Department of Commerce US Departmenl of Defense 
Change #4, January 1 1974 
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FIGURE 44 SAMPLE 
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FIGURE 45 SAMPLE FO‘ M MF3-31 B 
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FIGURE 46 SAMPLE FORM MF3-31 C 
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FIGURE 47 EXAMPLE OF COMPUTERIZED RADIOSONDE DATA 
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summarized in a teletype radiosonde observation (RAOB) message This message contains the 
meteorological and wind data at standard pressure levels (surface, 1000, 850, 700, 500, 400, 
300, 200, 150, 100, 70, 50, 30, 20, 10, 7, 5, 3, 2 and 1 millibar) and other levels as required by 
a particular mission These data are routinely sent to the World Data Center at Asheville, N C , 
for the SAGE project, however, copies of appropriate radiosonde data will also be sent to the 
Wallops Flight Center ASRP for distribution and archiving Typically, these data should be 
available at WFC within ten calendar days after acquisition for further handling and input into 
final data packages 

Standard card and plotting formats for comparisons of SAGE and correlative profile data 
are described in Section 4 2 5 


4.2. 4. 2. Super Loki Datasonde 

The Super Loki datasonde" is a rocket payload launched in support of the rocket ozone 
systems, primarily to provide temperature, pressure, and wind profiles The Datasonde payload 
is ejected at apogee to descend on a decelerator Temperature measurements are made between 
70 and 20 km during the descent phase Figure 48 shows the data flow at Wallops for this sys- 
tem The payload is tracked by an AN/GMD for the temperature/pressure data and by radar 
for positional data Both tracking systems at Wallops will digitize their data during the track 
The data flow at sites away from Wallops may be slightly different in that the AN/GMD data 
may not be digitized but may be a TMQ-5 paper strip chart The radar data may also be in ana- 
log form 

The digitized AN/GMD and radar data from Wallops, along with datasonde calibration 
data and lower-level balloon radiosonde data are inputted to the METROC-K computer program 
for data reduction Atmospheric density is computed by using the temperature data, a base-level 
pressure, the hydrostatic equation, the equation of state, and a standard baseline temperature 
Atmospheric pressure is also computed from this information, and winds are derived from the 
radar data Figure 49 is a typical set of datasonde data Figure 50 (the portion above the 26-28 
km level) is a typical chart presentation of the datasonde temperature/altitude data That por- 
tion below 26 km is radiosonde data 

For those cases in which the AN/GMD and/or radar data from other sites arrive at Wal- 
lops in analog form, a different approach is used The TMQ-5 meteorological data will be 
extracted manually and put through a FORTRAN formatting process The radar data will be 
digitized and then inputted together with the meteorological data to the Hypso-2 computer pro- 
gram at Wallops 

The meteorological data products from either the METROC-K or Hypso-2 programs are 
provided (as shown in Figure 51) as part of the supporting data required for processing the 
ozone data from the Super Loki and Super Areas ozonesondes 

A standard card and plotting format for comparisons of SAGE and correlative data is 
described in Section 4 2 5 


' ^Federal Meteorological Handbook No 10, Meteorological Rocket Observations, NASA, U S Department of Commerce, U S Depart- 
ment of Defense Change#! May I, 1977 
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4 2.4 3. Upper-Air Synoptic Analysis for 5-, 2-, 1- and 0.4-Millibar Surfaces 

Meteorological rocketsonde and satellite radiance data are used to generate high-altitude 
synoptic charts Broad-scale analyses for the Northern Hemisphere at 5-, 2-, 1-, and 0 4-mb lev- 
els are prepared by NOAA routinely on a weekly basis The WFC ASRP Office will coordinate 
the specific requirements for these charts with NOAA 


4 2 4 4. NMC Meteorological Data Products Provided Regularly for SAGE 

In addition to the above data provided by specific ground truth experiments, the Upper 
Air Branch (UAB) of the National Meteorological Center (NMC) will provide for each SAGE 
tangent location and time meteorological data profiles interpolated from NMC gridded global 
data sets, as well as profiles at the nearest radiosonde location Tapes containing these data will 
be sent regularly by NMC to NASA Langley Research Center for use in SAGE data validation 
or in routine data processing 

The starting point for deriving the profiles at SAGE tangent locations and times is an 
archive tape of gridded global height and temperature fields at constant pressure levels Each 
file on this tape contains one calendar week of fields in the format specified by NMC Office 
Note 84 Although these gridded archive fields are created weekly, they will be distributed on 
multifile tape in 2- to 4-week batches Figure 52 is an example of the log that accompanies 
each weekly archive file Note that 80 fields are currently available, in the vertical sequence 
shown in Figure 52, for each of the seven days The log serves only as a guide, since each 
available tape field can be identified through its unique 12-word label A more detailed descrip- 
tion of tape and field formats is available from NMC on request Interpolations within these 
fields are the source of all density-height profiles at SAGE tangent locations and times 

Figure 53 shows the four analysis programs that generate fields saved on the UAB archive 
tapes Input parameters to these various systems of analysis are also shown Within the next 
several months the transition from NOAA-5 (VTPR) satellite radiance data to TIROS-N 
fTOVS) data will affect all products at all levels A significant amount of work remains to be 
done in integrating TOVS data into the analysis schemes and profile error estimates 

The links between SAGE events, analyzed fields, density/height profiles, error estimates, 
and radiosonde data are diagrammed in Figure 54 The entire system of Figure 54 must run in 
quasi-real time (within ten days of the earliest SAGE event to be profiled in a given batch), 
since radiosonde data are available only for a limited time The profile-generation program will 
be run once each week (between Sunday and Wednesday) and will process batches of input 
information from the previous calendar week Figure 54 also shows input functions and pro- 
vides a general summary of the content of the output profile information, which will also be 
distributed in weekly batches 

Figures 55 and 56 and Table 13 depict specific processing of SAGE events, error esti- 
mates, and profile information, respectively Figure 55 indicates that an input tape of SAGE 
event cards will be accepted by UAB covering any length of time from several hours to several 
weeks, although the latter is preferred Each batch of SAGE events, as well as each event 
within a batch, is expected to be in chronological sequence The format of UAB’s analysis and 
data archives dictates the creation of profile information in calendar-week batches As in the 
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*Note that final format will have four additional lines (for a total of 80) Northern Hemisphere tropopause 
temperature and pressure, and Southern Hemisphere tropopause temperature and pressure 


FIGURE 52 EXAMPLE OF ARCHIVE LOG FOR TAPE FILE OF NMC GRIDDED GLOBAL DATA 
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FIGURE 53 NMC ANALYSIS ARCHIVE SOURCES AND THEIR INPUT 












FIGURE 54 


INPUTS AND OUTPUTS OF NMC PROFILE-GENERATING PROGRAM 
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FIGURE 55 EXAMPLE OF NMC PROCESSING FOR SAGE EVENTS 
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Table 13 


TENTATIVE FORMAT FOR NMC PROFILE OUTPUT 

Fourteen 80-character records are written for each event - all values 
can be read as integers or characters 


Parameter 

Units 

Fortran 

Format 

Remarks/Sample Element 

Event 

- 

(80A1) 

Original event card 

1000- 10MB Temp 

Tenths °K 

(1415, 10X) 

|a2732| = 273 2°K ^ 0°C 

5-0 4MB Temp 

Tenths °K 

(415, 60X) 

|a2732| = 273 2 o K*=0°C 

1000-0 4MB T ERR 

Tenths °K 

(18I4.8X) 

|aa56| = ±5 6°C error 

1000- 10MB HGT 

Meters 

(1415, 10X) 

|a 15801 = 500MB geopotential height 

5-0 4MB HGT 

Meters 

(4I5.60X) 

|42412| = 2MB geopotential height 

1000-0 4MBH err 

Meters 

(18I4.8X) 

|9999| = missing 

1000- 10MB Density 

gm/m 3 

(14(14, 11), 10X) 

15369(21 = 5369 X 10' 2 = 53 69 gm/m 3 

5-0 4MB Density 

gm/m 3 

(4(14,1 1),60X) 

|5112|4| = 5112X10' 4 = 5112 gm/m 3 

1000-0 4MB % Density ERR 

- 

(18I4.8X) 

|aa30| = 3 0% 

Radiosonde Ident, Tropopause 


(517,315) 

If available Stn, Date, Time, Lat, Lon, 
T, H, P 

1000- 10MB T r/s 

Tenths °K 

(1415, 10X) 

See above 

1000- 10MB H r/s 

Meters 

(1415, 10X) 

See above 

Miscellaneous 

- 

- 

- 


case of the gridded archive field files, profile information files will be distributed on multifile 
tapes containing several weeks of profiles 

Figure 56 gives an abbreviated depiction of the current estimates of temperature profile 
errors (e„) at mandatory pressure levels, assuming that events occur at 12Z Additional errors 
result from time interpolation tor event profiles between 12Z analysis cycles The equation 
e = e„(l + a |sin 7 t 7724|) supplies this adjustment error where a is 1 from 500-100 mb, a is a 
function of season from 70 mb through 10 mb, and a is 0 above 10 mb This scheme is 
presented in greater detail in Appendix D For the present, UAB is assuming that the relative 
error in density is nearly equal to the relative error in temperature In most cases the pressure 
error contribution to the density error is smaller than the temperature error contribution as 
defined above If the pressure error effect is determined to be significant, UAB may convert 
this error to an additional probable temperature profile error A slot is being allotted for report- 
ing height errors if they are later determined to be useful 
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A tentative format has been chosen for the SAGE event profiles and accompanying infor- 
mation Table 13 describes the proposed content of fourteen 80-character records that will be 
produced for each SAGE event (Note that more than 14 records may be used in the final for- 
mat ) This record count per event will remain constant, and any missing information will be 
filled in with an appropriate number of 9s Ample room has been left for revisions and addi- 
tions to the profile information Once format and content are agreed upon, such changes will 
be kept to a minimum to allow for a consistent data set 


4.3. Standard Card and Plotting Formats for Data Exchange 

To facilitate comparison the following standard formats are recommended for correlative 
profile data 


4.3.1. Punched Cards 

In general, punched cards are preferred over magnetic tapes as a data exchange medium, 
because of the possible problems entailed in reading tapes on one computer that have been 
written on another Because of the limited number of correlative profiles to be measured by 
any single investigator, the bulk of cards necessary for data exchange is not expected to be 
excessive Two types of card format are acceptable one with a single height, but with many 
parameters, per card and one with many heights, but with a single parameter, per card These 
alternative formats are described in Sections 4 3 12 and 4 3 13 Cards in either format must 
be preceded by main-deck header cards, as described in Section 4 3 11 


4. 3. 1.1. Main-Deck Header Cards 

The first two main-deck header cards are in (nAx) format and contain the following infor- 
mation 

1 Instrument name (cols 1-20), location name (21-40), date and Greenwich Mean Time 
(41-60), investigator’s serial number, etc (61-80) 

2 Beginning latitude and longitude (1-20), end latitude and longitude (21-40), investigator 
(41-60), institution (61-80) 

The third main-deck header card is in (15) format and gives the number of altitudes for which 
data will follow 


4. 3. 1.2. Data Cards: Single Height and Many Parameters per Card 

This type of format was used in the first SAM II ground truth experiment in Sondres- 
trom, Greenland, November 1978 It lends itself readily to cases with a limited number of 
measured parameters that can be anticipated sufficiently in advance 
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For aerosol data the standard FORTRAN coding format is 


(15, F6 2, F6 1, F5 1, F4 1, 1 PE9 3, OPF6 3*U*F4 3,F6 3*U*F4 3,F6 3*U*F4 3, 2F6 3) 
The allocation of fields to parameters is as follows 


Columns 

Formal 

Parameter 

Typical Sensor 

1-5 

15 

Card N umber 


6-11 

T6 2 

Altitude (km) 

Many 

12-17 

F6 1 

Pressure (mb) 

Balloonsonde, rocketsonde, aircraft 

18-22 

F5 1 

Temperature (C ) 

Balloonsonde,, rocketsonde aircraft 

23-26 

F4 1 

Dew-Point Depression (C) 

Balloonsonde rocketsonde, aircraft 

27-35 

IPF9 3 

Neutral Oas Density (g/cm 3 ) 

Balloonsonde, rocketsonde, aircraft 

36-41 

OPF6 3 

180° Scattering Ratio* 

Lidar 

42 

*U* 

Delimiter 


43-46 

F4 3 

Absolute Uncertainty (1 cr) in Preceding Parameter 


47-52 

F6 3 

25° Scattering Ratio* 

Polar nephelometer 

53 

*U* 

Delimiter 


54-57 

F4 3 

Absolute Uncertainty ( 1 <r) in Preceding Parameter 


58-63 

F6 3 

F.xlinction Ratio* 

SAGE, SAM II balloon photometer 

64 

*U* 

Delimiter 


65-68 

F4 3 

Absolute Uncertainty (1 cr) in Preceding Parameter 


69-74 

F6 3 

Number Mixing Ratio, r > 0 15 gm (mg -1 ) 

Dustsonde 

75-80 

F6 3 

Number Mixing Ratio, r ^ 0 25 gim (mg -1 ) 

Dustsonde 


At wavelength specified in header cards 


Fields for parameters not measured in a particular correlative experiment will be left blank 
The rationale for specifying ratio data (rather than absolute concentrations) is that this avoids 
using an E format and saves space on the cards By using the molecular density provided on the 
cards, the computer can readily convert to absolute concentiations 

For ozone and nitrogen dioxide data the standard FORTRAN coding format is 


(15, F6 2, F6 1 , F5 1 , F4 1 , 1 PE9 3, 0PF6 3*R*F4 2, F6 3*R*F4 2) 
The allocation of fields to parameters is as follows 
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Columns 

Format 

Parameter 

Typical Sensor 

1-35 


Same as for aerosol cards 

36-41 

OPF6 3 

0 3 Mass Mixing Ratio Gg/g) 

Balloonsonde, rocketsonde, photometer 

42 

*R* 

Delimiter 


43-46 

F4 2 

Relative Uncertainty (lcr) in Preceding Parameter 


47-52 

F6 3 

NO 2 Mass Mixing Ratio (m/ig/g) 

interferometer, spectrometer 

53 

*R' 

Delimiter 


54-57 

F4 2 

Relative Uncertainty (1 <x) in Preceding Parameter 



We anticipate correlative data that will not fit into the above formats For example, aero- 
sol number data may have different size cutoffs, and extinction or scattering data may be meas- 
ured at more than the allotted number of wavelengths For data that do not fit into the above 
format the following more general format is recommended 


4. 3. 1.3. Data Cards: Many Heights and Single Parameter per Card 

In this format each parameter has its own header card and data subdeck The complete 
set of subdecks follows the three header cards described in Section 4 2 5 1 1 The parameter 
header cards and data subdecks are as follows 

A Ititudes 

Altitude header card Format (nAx) Example 
ALTITUDE (MSL) — KM 
Altitude data cards Format (8F10 x) 

First Parameter 

First parameter header card Format (nAx) Example 

PARTICLE NUMBER GT 0 3 UM RADIUS — CM**(-3) 

First parameter data cards Format (8E10 x) or (8F10 x) 

First Parameter Uncertainty (tj available) 

First parameter uncertainty header card Format (nAx) Example 

ABS UNC , PARTICLE NUMBER GT 0 3 UM RADIUS — CM**(-3) 

First parameter uncertainty data cards Format (8E10 x) or (8F10 x) 

Header cards and data subdecks for other parameters measured at the same altitudes fol- 
low the first-parameter subdeck in similar format Parameters measured at different altitudes 
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must be preceded by three new header cards (as described in Section 4 2 5 1 1 ), plus a new 
altitude header card and altitude data subdeck 

Note that all height and parameter data described in this section can be read by the single 
FORMAT specification (8F10 x) or by the single specification (8E10 x) Thus, a plotting pro- 
gram with either of these input formats can read any data prepared in the format described in 
this section Such a program will be available at NASA Langley Research Center for plotting 
SAGE and correlative results on matching axes for rapid comparisons (See Section 4 4) 


4.3. 1.4, Other Formats 

It is to be expected that some correlative data will not fit into any of the formats described 
above (One example would be repeated measurements made at one altitude, but in varying 
horizontal positions) If this is the case, an appropriate format can be chosen by the correlative 
sensor scientist Cards in such nonstandard formats should be accompanied by FORMAT cards 
suitable for reading the data Brief instructions or a sample FORTRAN card-reading program 
would also be desirable 


4.3.2. Plotting Axes 

The following standard plotting axes are recommended semilogarithmic, 7 inch (17 78 
cm) x 10 inch (25 4 cm), four orders of magnitude (logarithmic) on the 10-inch dimension, 
and 7 major divisions on the 7-inch dimension The logarithmic axis is for constituent concen- 
tration or mixing ratio data, and the linear axis is for height data. The standard height ranges 
are 


Ozone Data 10-45 km Nitrogen Dioxide Data 15-50 km 

Aerosol Data 5-40 km Total Density Data 10-45 km 

Any limits and units can be used on the logarithmic scale Comparisons will be made by over- 
lapping plots and by sliding the logarithmic scales to obtain overlap 

As regards the plots, it is preferred that constituent data be expressed in absolute units 
(not mixing ratios), plots of mixing ratio, however, in addition to absolute concentration plots, 
will be welcome and useful 

Correlative sensor scientists (especially those outside of the US) who are unable to 
obtain the proper graph paper can request it from Mr Leonard McMaster, the SAGE Science 
Manager, at Mail Stop 234, NASA Langley Research Center, Hampton, Virginia 23365 


4.4. Joint Evaluations and Visits to Langley Research Center 

The SAGE Experiment Team (SET) feels that cooperative efforts by correlative and 
SAGE scientists are required to obtain the best scientific judgment of data validity This 
cooperative evaluation must be based on a thorough comprehension of the error sources in the 
correlative measurements Hence, a writeup explaining these sources and their effects should 
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be furnished the SET as part of the joint evaluation. The most useful way of expressing these 
effects is to include error bars on the correlative data points 

Communication would also be facilitated by having correlative and SAGE sensor scientists 
work together at a common location Toward this end, NASA plans to make space and faclities 
available at the Langley Research Center (LRC) for visiting scientists to analyze and evaluate 
SAGE and correlative sensor data Correlative sensor scientists are encouraged to visit LRC 
and take advantage of this opportunity to exchange information (Any such visits must be 
coordinated with NASA International Affairs Division, Code 21C-17, Washington, DC 
20546) 
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Appendix A 


THE SAGE INSTRUMENT AND EXPECTED MEASUREMENT ERRORS 


The SAGE sensor is a four-spectral-channel radiometer that measures the extinction of 
solar radiation during solar occultation As the conveying spacecraft emerges from the earth’s 
shadow during each orbit, the sensor will acquire the sun and measure solar intensity in four 
wavelength bands centered at 0 385 ftm, 0 45 0 60 /nm, and 1 0 /itm As the spacecraft 

continues in orbit, the line of sight from the spacecraft to the rising sun will scan the earth’s 
atmosphere, resulting in a measurement of the attenuated solar intensity at different atmos- 
pheric layers The procedure will then be repeated in a reverse sense during spacecraft sunset 

Each sunrise and sunset event will be monitored from the top of the clouds to approxi- 
mately 150 km above the earth’s surface The sensor will have an instantaneous field of view 
of approximately 0 5-arc minutes in elevation, which corresponds to approximately 0 5 km 
measured at the horizon for a 600 km orbit The total field of view is 360° in azimuth and -14° 
to -30° in elevation The dynamic range of each radiometric channel is approximately 3000 and 
the uncertainty in any radiometric measurement is specified to be less than 0 1% of the unat- 
tenuated solar intensity The sensor is partially self-calibrating, as a measurement of the unat- 
tenuated solar intensity is made prior to each spacecraft sunset and following each spacecraft 
sunrise Figure A-2 shows the orbit and viewing geometry of the SAGE instrument in orbit 

The instrument module (Figures A-2 and A-3) consists of optical and sensor subassem- 
blies mounted side by side The optical subassembly consists of a flat scanning mirror, 
Cassegrain optics, and a detector package The entire optical subassembly is gimballed in 
azimuth The azimuth servo employs sun sensors driven to null on the center of the sun to a 
tolerance of ± 45-arc seconds At the beginning of a sunrise or sunset event, the instrument 
slews in azimuth to a position for acquiring the sun Upon acquisition in azimuth the mirror 
servo scans in elevation until the sun is acquired The scan range is then constrained to scan 
back and forth across the solar image only 

The solar input is reflected from the scan mirror through the Cassegrain telescope, which 
produces a solar image at the science detector aperture This image is scanned across the aper- 
ture by the motion of the scan mirror The radiation through the aperture is dispersed by a 
holographic grating and the four beams representing the four wavelength bands are then col- 
lected and applied to silicon PIN diode detectors The outputs of the detectors are fed to 
signal-conditioning amplifiers whose outputs go to the PCM encoder 

After multiplexing and digitizing the signals, the PCM encoder transfers the digital data to 
the AEM-B data system The radiometric data for each wavelength channel will be sampled 64 
times per second or approximately four times per km of tangent altitude, and digitized to 12 
bits (3072 bps), these data, plus science-supporting data and instrument module housekeeping 
data, total 5440 bps 
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SURVEY OF POTENTIAL GROUND TRUTH SUPPLIERS 
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Appendix B 


SURVEY OF POTENTIAL GROUND TRUTH SUPPLIERS 


In February 1978 a letter and questionnaire (reproduced on the following pages) were 
composed to aid in locating established teams that might be able to provide correlative meas- 
urements for SAGE data validation Responses describing potentially useful measurements are 
summarized in Table B-l Many of the teams listed are now being integrated into the Euro- 
pean, Japanese, or United States SAGE ground truth groups As an aid in the scheduling of 
correlative measurements, all teams on this list will be provided with updated information on 
SAGE coverage and the activities of the established ground truth groups 
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The Stratospheric Aerosol and Gas Experiment (SAGE) is scheduled 
for launch in February 1979 aboard the National Aeronautics and Space 
Administration's AEM-B satellite. The purpose of SAGE is to map verti- 
cal profiles of ozone, aerosol, nitrogen dioxide, and Rayleigh molecular 
extinction. We expect that the ozone data will extend from about 10-45 
km, the aerosol data from cloud tops to about 35 km plus occasional 
strong layers in the mesosphere, the nitrogen dioxide data from about 
25-40 km; and the Rayleigh molecular extinction from about 15-40 km. 

SAGE is a 4-channel photometer that measures the intensity of sun- 
light (centered at wavelengths 0.385, 0.45, 0.60 and 1.0 pm) traversing 
the earth's limb during spacecraft sunrise and sunset events. In this 
manner it will measure vertical profiles of 4-wavelength extinction, at 
the rate of about 30 profiles per day Spatial coverage will extend 

from about 75° N to 75° S latitude (with some seasonal dependence) and 
thus will complement the coverage (64° - 80° N and S) of the SAM-II 
stratospheric aerosol sensor to be flown on the Nimbus G Satellite. 

SAGE's four-channel extinction measurements will be numerically 
inverted to yield vertical profiles of aerosol extinction (and inferred 
number density), ozone concentration, nitrogen dioxide concentration, 
and total molecular density (When available, molecular density may be 
derived from the rawinsonde network and other sources and used as an 
input to the 4-channel inversion process, if this improves the ac- 
curacy of the other derived parameters ) The derived data will be 
archived and made available to the scientific community for use in a 
variety of studies 


A very important 
tific community is to 
SAGE measurements and 
density measurements 
a primary function of 


step in making these data available to the scien- 
validate them by making comparisons between the 
other aerosol, ozone, nitrogen dioxide, and 
made nearby in space and time This validation is 


the SAGE Experiment Team 


M P McCormick (Team 


SRI International 

333 Ravenswood Ave • Menlo Park, California 94025 • (415) 326-6200 • Cable STANRES, Menlo Park • TWX 910-373-1246 
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Leader), R. A. Craig, D. E. Cunnold, G. W. Grams, B M Herman, D E 
Miller, D. G. Murcray, T. J. Pepin, W. G. Planet, and P. B. Russell . 

One of our first actions is to identify those data available from other 
measurements that can be used for validation purposes. Not only aerosol, 
ozone, nitrogen dioxide, and density measurements per se, but also 
related information, such as volcanic eruption and noctilucent cloud 
observations, will be useful to us. 


We would appreciate your assistance by filling out the enclosed 
questionnaire and returning it at your earliest convenience. Attached 
to the questionnaire are guidelines for the type of information required 
for our planning. Even if you do not plan to make correlative measure- 
ments during the SAGE validation period, please answer Part I regarding 
possible data use , as this will aid us in our efforts to accelerate data 
utilization. If, by returning Part II, you indicate an intention to 
make correlative measurements, we will contact you to discuss SAGE 
intercept times and dates for your sensor location(s). 

We look forward to your support of this very important phase of a 
major stratospheric research endeavor. 


Sincerely, 


Richard A. Craig 
Department of Meteorology 
Florida State University 
Tallahassee, FL 32306 


Walter G. Planet, S321B 
National Environmental 
Satellite Service 
NOAA 

Washington, D.C. 20233 


Philip B. Russell 
Atmospheric Sciences 
Laboratory 
SRI International 
Menlo Park, CA 94025 


(SAGE DATA VALIDATION GROUP) 
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PART II: AVAILABILITY OF DATA TO VALIDATE SAGE MEASUREMENTS 


The following measurements or observations, relevant to SAGE validation, 
are contemplated: (See sample responses on following cage) 

1. Parameter to be measured: 

2. Accuracy of measurement: 

3. Altitude region: 

4. Altitude resolution: 

5. Observation period: 

6. Frequency of observation: 

7. Measurement technique: 

8. Instrument type: 

9. Measurement platform: 

10. Type of data product: 

11. Funding authority: 

12. Measurement program: 

13. Status of prime instrument: 

14. Assurance of instrument availability m Jan. 1979-July 19 7 9 * 

15. . Location of measurement 1 

16. Experimental limitations: 

17. Instrumental physical characteristics: 

18. Experimenter operation experience: 


Dace: 


Return to : 


Prepared by Dr. Philip 3. Russell, 12056 

Atmospneric Sciences Laboratorv 

Adcress: SRI International 

Menlo Park, CA 94025 
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Guidelines for Parc II 


The following guidelines are co assist you in supplying che 
necessary information for Parc II of che questionnaire. 


Required Parameter 

1. Parameter to be measured: 

2. Accuracy of measurement: 

3. Altitude region: 

I. Altitude resolution: 

3. Observation period: 

5. Frequency of observation :- 

7. Measurement technique: 

8. Instrument type: 

9. Measurement platform: 

10. Type of data product: 

11. Funding authority: 

12. Measurement program: 


13. Status of prime instrument: 

11. Assurance of instrument of 

availability, Jan 1979-July 1979 
13. Location of measurement: 

15. Experimental limitations: 

17. Instrumental pnysicai 
characteristics 
13. Experimenter operation 
experience: 


Example of Responses 

Aerosol particle number: Ozone 
concentration; noctilucent cloud 
occurrence 

9f , — 1 — ] 

/• , ± a m sr x 
a km to b km 
c km 

June-August 1979 
single flight weekly 
in situ, remote 
lidar, sampler 
RB-57 at 20 km altitude 
direct, analysis 
FAA, COVOS 

Atmospheric radiation studies, 
part of WWW, flight test for 
specific instrument development, 
etc . 

x years of use, demonstration 
planned for : rasearcn stage 

” prooability 

Siberia; Fairbanks, Alaska 
day-night, duration 

weight, size, power requirements 

10 flights m R3-37 
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Table B-l 


CORRELATIVE MEASUREMENT CAPABILITIES DESCRIBED 
BY QUESTIONNAIRE RESPONSES 


Part A Ground-Based 


Team 

Objective 

Instrument 

Site 

CSIRO,Aspendale 
(R Kulkarni) 

1 Total 0} 

2 O 3 profiles 

1 Dobson spectrophotometer 

2 Dobson Umkehr observations 

Aspendale 

Vairns 

Brisbane 

Perth 

Hobart 

Macquarie Island 

NCAR 
(C Frush) 

Aerosol profiles 

Ruby lidar* 

(X = 0 6943 /im) 

Boulder, 

Colorado 

INPE, Sao Paulo 
(B Clemesha) 

Aerosol profiles 

Dye lidar 
(X= 0 5893 nm) 

Sao Paulo, 
Brazil 

(23 S, 46 W) 

U of West Indies 
(G Kent) 

Aerosol profiles 

Ruby lidar 
(X= 0 6934 fim) 

Kingston 

Jamaica 

Hebrew U of Israel 
(A Cohen) 

Aerosol profiles 

Ruby lidar 

Jerusalem, 

Israel 

Atmos Env Serv 
(J Kerr) 

Total NO; 
above 14 km 

Twilight 

spectrophotometer 

Downsview, 

Ontario 

(44 N, 79 W) 

U of Arizona, 
(R Schotland) 

O 3 profiles 

Dye lidar 

Tucson, AZ 
(32 N, 110 W) 


Part B Airplane-Borne 


Team 

Objective 

Instrument 

Site 

NASA Ames 
(M Lowenstein) 

0 3 ,N0; 

concentration 

Chemiluminescent 
sensor on U-2 or 
other aircraft* 

U-2 routes, based at 
Mountain View, CA 
(37 N, 122 W) 

NASA Ames 
(N Farlow, 

G Ferry) 

Aerosol size 
distribution, 
concentration, 
composition 

Impactor on U-2,* 
Lear jet, or balloon 

California, Alaska 
Panama, Hawaii 

LASL 

(W Sedlacek) 

Aerosol concentration 
and composition, 

4 altitudes 12-20 km 

Filter sampler and 
A N counter on 
WB-57F 

0-75 N 

Panama-Houston- 

Seattle-Alaska 

NASA Lewis 
(D Rriehl) 

Aerosol concentration 
and size distribution, 
6-14 km 

Forward-scattering 
counter on commercial 
B-747’s (GASP) 

United and Pan/ Am 
airline routes 


‘Already included in the U S ground truth experiment schedule See Table 2 
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Table B-l (Concluded) 


Part C Balloon-Borne 


Team 

Objective 

Instrument 

Site 

U ofKoln* 
(A Ghazi) 

0 3 , HjO profiles, 
0-35 km 

Optical sonde 

Koln, W Germany 

CSIRO, Aspendale 
(R Kulkarm) 

O 3 profiles 

Mast sonde 

Aspendale 

CSIRO, Aspendale 
(I Galbally, 

with Canadian AES and 
York U ) 

O 3 ,NOj , NO profiles, 
10-35 km 

Chemiluminescent 
and infrared 
spectrometer 

Mildura (34 S) 

Alice Springs (23 S) 

U of Wyoming^ 
(J Rosen) 

Aerosol, 
H,,0,03,CN 
profiles, 0-30 km 

Optical particle 
counter 
(plus others) 

Laramie, WY 
(41 N, 105 W) 
Arctic, Antarctic, 
Equatorial 

AES, Ontario 
(W Evans) 

O 3 ,NOi ,NO 
HNO 3 , Aerosol 
H 3 O profiles 
10-36 km 

Stratoprobe 
Payload (Several 
instruments assembled 
for LIMS/Nimbus G 
rendezvous) 

Cold Lake, Alberta 
(55 N, 110 W) 
February 1979* 

JPL 

(C Farmer) 

O 3 ,N0 2 ,other 
profiles, 

tropopause-42 km 

Fourier 

spectrometer 

Palestine, Texas, 
Alice Springs, 
Australia 


‘This team has recently joined the SAGE European ad hoc ground truth group 
^Already included in U S ground truth experiment schedule See Table 2 
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Appendix C 

INFORMATION ON NOCTILUCENT CLOUD OBSERVATIONS 
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Table Cl 


DESCRIPTION OF GROUND-BASED PASSIVE NOCTILUCENT CLOUD OBSERVATIONS 

1. Parameter to be measured: Noctilucent cloud occurrences; time of display, 

intensity; forms; extent of azimuth; extent in elevation. 

2. Accuracy of measurement: Several degrees in azimuth and elevation. 

3. Altitude region: 'V 83 ± 10 km 

4. Altitude resolution: Not actually measured, but see attached Fig. Cl 

from "Noctilucent Clouds" by V.A. Bronshten and N.I. Grishin, Keter Publ., 
Jerusalem. 

5. Observation period: March 1 - October 31, annually. 

6. Frequency of observation: Daily 

7. Measurement technique: Visual estimation (supplemented by theodolites at some 

locations); photographic recording of displays at 6 stations. 

8. Instrument type: Human eye; 35 mm camera 

9. Measurement platform: Surface of the earth. 

10. Type of data product: Direct, recorded on forms for punching onto cards; 
printed in annual publication 

11. Funding authority: Atmospheric Environment Service (Canada) 

12. Measurement program: International NLC program in cooperation with World 

Meteorological Organization 

13. Status of prime instrument: Operational for many years 

14. Assurance of instrument availability in Aug. 1978-Sept. 1979: 1007. 

15. Location of measurement: 60 stations in Canada; 16, in U.S.A. 

16. Experimental limitations: Observations made during pre-sunrise and post- 

sunset periods when sun's depressional angle is between 6° and 18°. 

17. Instrumental physical characteristics: Normal as to class 

18. Experimenter operation experience: Annual regular observations since 1964 

Prepared by: Mr. E. J. Truhlar 

Address. Atmospheric Environment Service 

4905 Dufferin Street 

Downsview, Ontario, Canada M3H 5T4 

Date: May 18, 1977 
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FIGURE C-1 



DISTRIBUTION OF NOCTILUCENT CLOUD HEIGHTS FROM 695 MEASURE- 
MENTS BETWEEN 1887 AND 1964 
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Your file Votre dossier 


Environment Environnement 
Canada Canada 


Atmospheric Environnement 
Environment atmospherique 

4905 Dufferin Street 
Downs view, Ontario 
M3H 5T4 


Our hie Notre dossier 8061-2 (ARPD) 


May 5, 1977 


Dr. P.8. Russell 

Science Coordinator for Groind Truth 
SAM-II Nimbus G Experiment Team 
Stanford Research Institute 
Menlo Park. Calif. 94025 
U.S.A. 

Dear Dr. Russell : 

This letter will confirm preliminary arrangements made during 
your recent discussion with Mr. E.J. Truhlar concerning the provision of 
noctilucent cloud (NLC) observations in support of SAM-II measurements of 
mesospheric aerosols. Observations of NLC data will be forwarded to you 
after they have been received from the stations, transferred to cards, pro- 
cessed by computer for quality and listed in tabular format. (See attach- 
ment for an example of a regular listing of such data, including the explan- 
atory legend). About a one-month delay should be expected before a listing 
for a particular data-month is received by your institute. 

The following table shows the average monthly distribution of 
station-night sightings of NLC during the period 1964-1976 inclusive. 

MAR APR MAY JUN JUL AUG SEP OCT 

0.6 0.8 2.6 40 71 24.5 1.2 0.6 

Most occurrences are in the months of June to August, with a peak in July; 
very fewoccur from March to May or in September or October. No NLC are observed 
from November to February during the fall and winter when the sun's elevation is 
too low to allow the clouds to be illuminated in the pre-sunrise and post-sunset 
twilight periods. 

We would appreciate receiving information on the progress of the arrange- 
ments to implement the SAM-II project. 


Yours sincerely 



B.W.Boville, Director 

Atmospheric Processes Research Branch 


145 



Table C2 
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(See next page for key to column headings.) 
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Table C2. 
(continued) 

KEY TO COLUMN HEADINGS 


- Year, month, night 

- Station 

- Station Identifier 

- Latitude 

- Longitude 

- Greenwich Mean Time 

- Local Solar Time (or Local Apparent Time) 

- Solar depression angle at time or sighting 

- NLC intensity on 5-point scale, from very weak to 
extremely bright 

- Structural forms : 

1 veils; 2-bands; 3-billows; H-whirls; 5-amofphous 

- Azimuthal extend of NLC, relative to geographic Nortn 

- Extent of NLC in elevation, relative to the horizon 

- Troposoheric cloud cover in twilight section of sky 
A-clear; 3-scattered; C-broken; D-overcast 
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HEAD OF DEPARTMENT 
D H MclNTOSH, BSc, M A , OSc , F R S E 

TEL 031-667 1081 Ext 2920 



Dr. P.B. Russell, 

Science Coordinator for Ground Truth, 
SAM- II Nimbus Experiment Team, 
Stanford Research Institute, 

Menlo Park, 

California 94024, 

U.S.A. 


DEPARTMENT OF METEOROLOGY 
THE UNIVERSITY 

JAMES CLERK MAXWELL BUILDINI 
KING'S BUILDINGS 

EDINBURGH EH9 3JZ 


7th October 1977 


Dear Dr. Russell, 

Thank you for your letter of 3rd October expressing interest 
in data relating to noctilucent clouds. 

We shall of course be pleased to cooperate with you in any way 
we can be sending you our data. Observations made here or reported 
to us are in large measure confined to June and Julyj only exception- 
ally are the clouds seen by ’our’ observers in August (or May). It 
seems likely then that our first reports of interest to you will be 
for 1979. 

The latitude belt you refer to (64° - 00°) is a good deal poleward 
of our most northerly observers. We shall make enquiries as to how 
far poleward we may be able to extend our network, perhaps obtaining 
the cooperation of other observers in Scandinavia or Iceland. 

Yours sincerely. 



/ (/ /jo • 




D.H. McIntosh. 
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Appendix D 

GROUND TRUTH PLANS FOR SAGE 
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Appendix D 

APPROACH FOR ESTIMATING ERRORS IN DENSITY PROFILES 


We describe here briefly the rationale for estimating errors in the temperature profiles to 
be required by SAGE and enumerate these errors for specials sets of conditions Constant- 
pressure synoptic maps will serve as the basis for the profiles The maps themselves are subject 
to several types of errors, viz, instrument errors, first-guess errors, and analysis-system errors 
It is known, in a general way, how these errors vary with geographical location, instrument 
type, altitude, and density of data It is also known approximately what errors are induced by 
interpolating between map times (errors due to moving weather systems between observations) 

The numerical values given here as error estimates result from long experience on the 
part of the Upper Air Branch staff in dealing with problems of instrument accuracy and data 
compatibility A complete rationale for these estimates is much too complicated for inclusion 
in a brief document, but this may be found in the papers cited among the references to this 
appendix 

The Upper Air Branch (UAB) of the National Meteorological Center (NMC) will provide 
atmospheric temperature profiles and associated error estimates for any point on the globe in 
support of SAGE operations All mandatory levels between 500 and 1 mb (~ 5 to — 48 km) 
are to be included in the profiles 

The principal systems of upper air observations for direct use in this support are rawin- 
sonde and satellite (Rocketsonde data comprise an essential contribution to the derivation of 
satellite profiles, however, besides providing "ground truth" for high-level radiosonde observa- 
tions ) Since spatial and temporal separations between the data points are large for all these sys- 
tems, it is necessary to use methods of interpolation, in both space and time, to derive the 
required temperature profiles The synoptic meteorological charts indicating temperature, wind 
and geopotential heights for constant-pressure surface areas serve as the basis for this interpola- 
tion 


The procedure for supplying temperature profiles may be outlined as follows 

(1) The SAGE team notifies UAB of its need for a temperature profile 
and temporal-geographical coordinates (r, $, X), where x = time, 
<f> = latitude, and X = longitude 

(2) The UAB consults all constant-pressure charts from 500 mb to 1 
mb valid for the map time Xj prior to the given time x, as well as 
the map times Tj subsequent to x The charts provide the to 
mates of temperatures for times X| and x ie, for (xl ’ 
and (x^, <t > , and X) Linear interpolation in time is then used to 
obtain the profiles for (x, d>, and X) Admittedly, there are errors 
associated with this profile, such errors are functions of time and 
pressure level 

(3) The UAB attaches "error bars" to the temperature profiles in accor- 
dance with the procedure outlined in the following paragraphs 
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When error estimates are made, the atmosphere between 500 mb and 1 mb is stratified in 
the vertical according to circulation regimes and measurement techniques The upper tropo- 
sphere, because of its dynamic features and the fact that it is monitored quite well by rawin- 
sondes, is taken here as a distinctive layer (500-100 mb) The lower stratosphere (50-10 mb) is 
the next highest layer, taken as distinct from the layer below because of its thermodynamic 
structure (Staff, Upper Air Branch, 1967, 1969) The upper stratosphere (5-1 mb) is the third 
layer, regaded as distinctive primarily because of the different techniques used for its measure- 
ment (Staff, Upper Air Branch, 1975) 

Figure D-l is a flow diagram that shows how errors in a given tropospheric temperature 
profile may be deduced for the Northern Hemisphere The first consideration is whether the 
number of rawinsonde or satellite data points within a distance a (radius of the scan circle A) 
exceeds a critical value x that depends on a If the answer is "yes," A is considered to be a 
data-rich area and, as the figure shows, the combination of radiosonde/satellite data determines 
the error characteristics— since any analysis system presently used would provide an analysis 
faithful to the mean values within a data-rich area Similarly, if the answer is "no,” A is con- 
sidered to be a data-poor region, and any errors would depend more on the analysis system 
used (including first-guess fields) than on current data 

The standard error of the temperature estimate at any point in space-time will thus 
depend on instrument accuracy and quality of analysis, as well as the time-interval t which 
separates the instant for which the profile is required from the time of the map on which the 
estimate is based The standard error e„ at time t = 0 is known from various studies of instru- 
ment performance and data compatibility (Mclnturff and Finger, 1968, Finger et al , 1973, 
Mclnturff, 1978) The empirical formula 


1 +a 


sin 


irt 

24 


where a = constant for a specified latitude, geographical area, and season, relates the sandard 
error at zero time to the standard error at an earlier or later time (the time-interval t is meas- 
ured in hours) (Bengtsson, 1969) The choice of a = 1 for the layer 500-100 mb is based on 
the observation that the error approximately doubles in 12 hours under the average conditions 
of variability that characterize the troposphere 

Figure D-2 is a flow diagram that shows how errors in Northern Hemisphere profiles 
between 70 mb and 10 mb may be deduced The situation for the Northern Hemisphere is 
simpler in this case, because VTPR data are not presently being used in the analysis (although 
satellite data, which will form the basis for such Southern Hemisphere analyses, may be used 
for these levels in the Northern Hemisphere as well) The standard errors vary seasonally more 
than in the troposphere, because the contrast between summer and winter is greater in the stra- 
tosphere As can be seen from the formulas, time interpolation should not pose much of a 
problem for the stratosphere in summer The wintertime stratosphere, however, is as variable 
as the wintertime troposphere, and this fact is reflected in assigning a — 1 for winter 

Even though there are in situ data available for the layer from 5 to 1 mb (primarily rock- 
etsonde data), they cannot be depended upon for present purposes as they are very sparse both 
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IS 

C(VTPR U R/S) 
> x 7 



YES NO 



e 0 = 1 2C 


e = e 


o 


( 1+ | S ' n ^l) 

KEY 


e = STANDARD ERROR AT TIME t 
= STANDARD ERROR AT t - 0 

t - DIFFERENCE BETWEEN TIME OF 
SATELLITE SOUNDING AND MAP 
TIME (in hours) 

R/S RADIOSONDE 

C(VTPR) COUNT OF VTPR REPORTS 
WITHIN SCAN REGION A 

C(R/S) COUNT OF RADIOSONDE 
POINTS WITHIN SCAN AREA A 

(VTPR U R/S) COLLECTION OF 
VTPR AND R/S POINTS 

CIVTPR U R/S) COUNT OF VTPR 

AND R/S POINTS WITHIN REGION A 

A CIRCULAR SCAN REGION OF 
REGION a 

P DATA-POOR REGION 

R DATA-RICH REGION 

x INTEGER DEPENDENT ON a 

A C R "REGION A IS INDLUDED 
IN REGION R" 

A C P "REGION A IS INCLUDED 
IN REGION P" 


FIGURE D-1 FLOW DIAGRAM OF ERROR-CALCULATION PROCEDURE FOR THE NORTHERN 
HEMISPHERE, 500- to 100-mb HEIGHTS 
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FIGURE D-2 



e„ = ± 1-2C e = ± 2-3C 

o o 


e = e Q (1 + alsin jrt/24!) 
a = 0 4 in summer (May-September) 
a = 0 5 in April and October 

a = 0 5 in winter (November-March) south of 30° N 
a - 1 0 in winter (November-March) north of 30° N 

FLOW DIAGRAM OF ERROR-CALCULATION PROCEDURE FOR THE NORTHERN 
HEMISPHERE, 70- to 10-mb HEIGHTS 
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in space and time Thus, satellite data derived with a currently used system will form the basis 
for producing temperature profiles Estimates of these errors are as follows 


At 5 mb, e ~ ±6C 
At 2 mb, e — ±8C 
At 1 mb, e — ±9C 

No time dependence is indicated— not because it does not exist, but rather because it is not 
known (It is believed that the error due to off-time data may commonly be as great with t = 1 
hour as with t * 10 hours, cf Miller and Schmidlm, 1971 ) 

As is well known, TIROS N is expected to supersede VTPR in the near future as the 
NO A A operational satellite For this reason we have included Figure D-3, which shows tem- 
perature standard errors in TIROS N simulations, both with and without noise It will be 
necessary to study the actual TIROS N data before Figures D-l and D-2 or the estimates of 
temperature errors given above can be modified to take the TIROS N contribution into account 

Most regions of the Southern Hemisphere (S H ) may be considered data-poor, except for 
the availability of satellite data The dearth of either S H rawinsonde or rocketsonde data 
makes calibration of satellite data difficult Consequently, error estimates for the S H are ques- 
tionable The following summarizes the best current error estimates for S H profiles 


Time dependence 


= e. 


1 +a 


sin 


El 

24 


where e 0 
t 

a 


error estimate for t = o 

interval separating map time from time for which profile 
is required (in hours) 

constant for a given altitude and geographic area 


Values of e 0 


Co 


e 0 


±2-3°C for 500-10 mb 
±6°C at 5 mb 
±8°C at 2 mb 
±9°C at 1 mb 


Values of a 

For 500-100 mb 

a =10 

For 70-10 mb 

a = 0 4 in summer (November-March) 

a — 0 5 in April and October 

a = 0 5 in winter (May-September), equatorward of 30°S 

a ~ 1 0 in winter, poleward of 30°S 

For 5-1 mb 

a =0 (time dependence unknown) 
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PRESS! 



STANDARD ERROR (TEMPERATURE) — °C 


FIGURE D-3 


TEMPERATURE STANDARD ERRORS IN TIROS-N SIMULATIONS 
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